The Design Criteria for 


Chrysanthemum Shells 


Takeo Shimizu 


Hosoya Fireworks Co. 


April 1955 


Table of Contents 


I. sPrefaées (NTs eet chances lccthe desecrate ee eateaeensed esaeant ee eadeceeat 
2. Construction of Chrysanthemum Shells oo... ccccccceseceseceseceeeceeecseeceeeeeeeeeeseeeseeeeeeaaes 
Problems in Designing Chrysanthemum Shells .........cccceccesceseceneceneceeeeeeeeeeeeeseeeseeenaes 
3.1. Factors that Have an Influence on the Initial Velocity of Stars oes 
3.2. Calculation of the Breaking Strength of a Shell oo... eeceeeereeeeceeeeeeeeeaeens 
3.3. The Effect of the Type of Paper on the Shell oo... cece ccecceesceeceeeeeeteeeteeeseeenaes 
3.4. Mass per Unit Cross-Sectional Area of a Spherical Star (Sectional Density) ..... 
3.5. The Compositions of the Burst Charge .........cccccccscsseceseceeseeeeceeeeeeeeeeseeeeeeseeeaaes 
3.6. The Force of Explosives f ......cccccccccsscsersesscssscesstecensceseceaesseccseccsesesesaseneeestesseeestes 
3.7. The Vivacity A of Burst Charge ......cccccccccccsscessecsseceseceseceeeeeeeseeeeeeeesseeeseeeeeeaaes 
3.8. Loading Density ccc cceccccccccssessseesseeeceeeecssecesecesecesecsseesseeseneeeneceseeeseeeeeeseeeaees 
3.9. Factors Concerning the Velocity of the Stars .......ccccesccesceeesceeeeeeeeeeseeeneeeteeenaes 
3.10. The Relation of Physical Conditions and the Visual Beauty 00... eeeeeeeeeeeees 
3.11. Calculating the Number of Stars in a Round Shell oo... eee eeeeseeseeeeceeeeneeeeeeeeens 


Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


Page 1 


Page 11 


Method for Experiments on Shell Burst, Preparation and 


Preliminary Calculations ........ccccccccscccsseesecessceseceeeceeecseceeeeeeeseeeseecseecsaecsaeceaeeeseenaeenseeeas 18 
4.1. Experimental Method oo... ccccecccssccssscessceeeceeeeeeeeesseeeseecseeesaecsaecsaecnseeesecneeneeees 18 
4.2. Aims of the Experiments ........cceccccccescessceseceeeceeeeeeseeeseecseeeseecsaecaecnsecnseceaeenseenas 18 
4.3. Types of the Sample Chrysanthemum Shells, their Construction and Preparation 18 


4.4. 


4.5. 


ASi1s Stars). sihactacaseindie eis eaves eel es ae eet 18 
(1) Types of Stars and their Symbols 0... ccceccseesseesteceteceteeeteeeseees 18 
(2) Construction Of Stars ....cccccceecceescesscesseeseesecsseceseceseceseeseeeeeeeeeneeaes 18 
(3) Composition of Stars ....ccccccccecsseesseeseeseesecssecsseceecseeseeeseeeeeneeens 19 
(4) Diameter of Stars oo... ccceccsseceseceseceseeeseeeseeceeceecseecssecesecneesseenes 19 

(5) The Mass, Sectional Density, and Density/Mass of one 
Titer OF StALS os ces cei yaugoencasies oeatenssecabsveagetencnsecubbeethgusedetann sates oan tek 21 

(6) The Burn Time of Stars in the Air when the Stars 

Are Not Moving ..ceccccesccesscesscssecseeceeeceeseesseecseeeeeeeecsecsaecnseenseenes 21 
A:3.2>- SBUrsUCHAarge. disscccestteiectscdenesiteeinbits dette dawnt iin tapecetes 23 
(1) Types and Symbols o....cececceecceescesseeeceeeeeseeeseeceeceecnsecsseenseenseenes 23 
(2) Construction of the Burst Charge .......cccecccecsseesseesteeeteceteeeteeneeees 23 
(3) Sensitivity test by Drop Hammer ou... ee eccccceesseesteesteceteeeteeeseeees 25 
(4) Calculation of Characteristic Values of Burst Charges ..............64 25 
1. Calorimetric Data of Solid and Gas oo... eeeeecceeseesteerteeneees 25 
a) The Internal Energy of Gases ........ccceccccsseesseeseeteeeteeeseeees 25 
b) The Internal Energy of Solids oo... cece eeceesseesseesteesteeneees 26 
2. The Characteristic Values of Bursting Charge (H) ............... 28 
The Characteristic Values of Bursting Charge (P)_ ............... 31 
. The Characteristic Values of Black Powder (S)_ ...........c0 33 
4.3.3. Fuse, Shells and Pasting o...c.ccccceccccsseescecscecsecsecseceseceseceeeseecneenaeens 35 
4.3.4. Construction of the Sample Shell oo... cece eccecsseesseereeseesteeseeeeenaeens 36 
4.3.5. Table of Sample Shells oo... cecccccccsscssccseeceeeeeeeeeeeeseecseeesecseceeeeaeens 37 
Preparation for the Test ......ccccccccccccsscssscssscsesceseeeeeeeeseesseeeseeeseecseecaecnaeenaecnaeenaeens 42 
4.4.1. Main Instruments .0.... eee ceeecesececeeseeseceeeeseeeecaecaeeeeceaecaeeeneesecneeeneeeees 42 
4.4.2. The Arrangement of the Instruments for Bursting Experiment .............. 43 
Progressing of the Experiment oo... cccccccccsseestecsteceeceeceeeeeeeeeeeeeeeeceseeeseeeaneeaaees 44 


Supplemental Experiments: Measurement of the Burn Velocity of 
Burst Charges and Stars in the Normal Atmosphere .........ccececesseesceseeeeeeeceseeeeeeeeaeenee 46 


Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


6. Formulae for Calculating the Velocity of Stars ......eeccceescescceseesceseeeceeeeaeeeeceaeeeeeneeaes 48 


6.1. Arrangement of Experimental Data 


6.1.1. Photography with a Revolving Shutter Camera .........ececeeeeseeseeeeeereeees 48 
6.1.2. Photography with a 16-mm Movie Camera .0....... ce eceecceseeeeerceteeneeeneeaee 48 
6.2. Trajectory of the Star in the Air oo... eee ccccsecsteceteceeceeeeeeeeeeeeeneeeseeeseeeeeesaes 48 
6.2.1. The Reaction of Air upon the Star oo... ceccesccesceseceseceeeeseeeeeeneeeseeesees 49 
6.2.2. The Burn Time of Stars oo... ecececsseecceesseeeeceseceeeeceeseceaeeaeeeeceaeeaeeereeaeeaee 70 
6.2.3. The Sectional Density of Stars oo... ccecccecsseesecesecseeeeeeceeeeeeseeeeeeseeeeeesaes 73 
6.2.4. Calculation of a Function F(v) = Av" oo. ecccecccccseescseesseseeseeeseseeeesseesenees 74 
6.2.5. Calculation of Values n and k for the Motion Equation 40. ou... 86 
6.2.6. The Influence of the Weather Conditions on the Ballistic Constants .... 88 
6.2.7. The Differential Equation Describing the Movement of Stars ............ 88 
6.2.8. Formula for a Star Moving in a Horizontal Direction 0.0.0... eee 88 
(1) For Stars that Consist of a Single Composition ..........ceceeeesteenees 89 
(2) For Stars that Consist of more than one 
Type of Composition (Color Changing Stars) ........ccceseeseeseeeees 90 
6.2.9. The Curvature of the Trajectory of the Star oo... ee eeceesseeteeteeeteeeees 92 
6.3. Imitial Velocity Of Stars oo... cecccsccesscseeceeseeeeeesseeeseecseecsaecaeceseceeeeeeeeeseseeeeseeeaaes 93 
6.3.1. The Formula of Initial Velocity Deduced from a 
Dot near the Bursting Point .........cccecccecsseeseesecsecneceseenseceseeeeeseeeseeeeses 93 
6.3.2. Calculation of the Velocities of Stars near the 
Bursting Point from the Dots on Photographs ..........cceccesseeeseeeteeeeeees 94 
6.3.3. Calculation of Initial Velocities at the 
Bursting Point from the nearest Dots of the Trajectories... cee 99 
6.4. Calculation of the Initial Velocity of Stars using a 16-mm Movie Camera 
and Comparison of those Results with those Previously Calculated 0.0... 105 
6.4.1. The Scale of the Photograph using the 16-mm Movie Camera ............. 106 
6.4.2. Number of Frames for Each Star Image ...........ccccescceseceseeeseeeseeeteeeseeesees 107 
6.4.3. Comparison of the Data for the Initial Velocities of Stars using 
the Revolving Shutter Camera and the Movie Camera ..........:esceeeeeeeee 108 
6.5. Initial Velocity of Star due to the Explosion (Bursting) ........c.ccccseesseeseetteenees 109 
6.5.1. The Influence of Various Manufacturing Processes on the 
Value of the Initial Velocity of the Star oo... ccc ecceesseeeseeseceeeeeeeeeeeeeees 109 
6.5.2. Comparison between the Initial Velocity for Ring Shells 
and: Full StarShell$:s..cicse.ssen pceasvertatecn sp taa pdetedssenterons wae oagiacteats 110 
6.5.3. Observation of the Burning of the Bursting Charge 0.0... eeeeeeteertees 111 


Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. Page ili 


6.5.4. The Equations for Calculating the Approximate Initial 


Velocity Of a: Stak: ssccasssassts seisiusishscasnasal desea steageactudeestesatusiekeateansenteet neds 114 
(1) Determination Of 06 ..ecccecccccccssecsseestecsteceteceseceseeseeseeeseeeseeeeeneeeaes 114 
(2) Determination of N (using Japanese Paper as the Standard) ......... 115 
(3) Determination Of Q ciecccccccccccccssecsseestecsseceteceeceeceeecseeeseeeseeeeeeeeaes 116 
(4) Determination of the Value Bo o..cecccccccccccseeseescsseeecseeseeseeseeseaees 118 
(5) Determination of Vivacity A for Each Bursting Powder ............... 119 
6.5.5. Adaptability of the Empirical Equation of Initial Velocity of a Star ..... 119 
(li) CalculatiomotK, cts icinctpptiethlteeasteixieall witteul arias aera ates 119 

(2) Determination of the Empirical Equation of the Initial Velocity 
OF A OUAL 24s, sateen edict vasatncss ceesneb ida shseats tes cet sa taveted seuss steateceternsout isan’ 124 
7. Equations for Practical Use in Designing Shells oo... cccccccesseesseceteceteceteeeteeeeeeeeneeees 126 
7.1. The Characteristics of Stars and the Preparation of Tables of Star Ballistics ..... 126 
7.2. Investigation of Star Ballistics in the Air from the Ballistic Tables... 126 
7.3. The Method for Obtaining the Function F\(1',4,,b) from F(7',A5,0)  ...eeceesceeereees 127 
7.4. An Examination of the Equations of Initial Velocity of Stars 0... eeeeeeeeee 128 

7.5. The Relationship between the Deviations of Initial Velocities AV/V and the 

Deviations of Flight Distances Axy/X ...cccccccccssceseceseceeeceeeeeeeeeeeeeeseeeseeeaeeesaeestees 129 
8. Examples of Designing Chrysanthemum Shells o.......ccceeceecesseeeneeereeeeeeteeeeeeseeeaees 131 
8.1. Comparison of Calculated and Experimental Values .........eceeeeceeseeseeteeneeereeee’ 131 
(a) Calculation of Initial Velocity oo... cc cecccsccesccessceseceeeceeeeeeeeeeseeeseeeeeeseeenseens 131 
(b) Calculation of the Flight Distance of the Star in the Horizontal Direction .. 132 
8.2. An Example of Designing a Chrysanthemum Shell... cee eeeeeeeeeeeneeeeeeeees 133 
9. Discussion and Conclusion .........esceescescessesseeseceseeseeeceaecaeeeeceaecaeeeceeseceeeseeseceaeeaeeereeaeeas 140 
References: ah btaSadave cas oSidonedbagathsdeadgtabedadardendsdieashieasnesduauacd odegassanagansavecdaedds 141 
Appendix © ebb thes rentcendatvancedeatodas cakte voebiyune seetatleadebganesv cues, vo libeerebeceubees eecae recede 143 
‘Tabled F5§(6= 0)  ~saeeetipieieaeiseedd eee vein eee eae eee eee 144 
Pable2rs DA tase Siti aS aceunn'gzs Sted aese An adatsmsheunegnesiavenl satdasangetivem, Mar cesetysaencynoss® 152 
‘Fable:3> fy (b= 0) cath deen Madea el te tee ont at a em aed es 153 
Table: (Fy (DS OS) ~cpecteschnshaGnceed ieee teh bociontee see ant eee se veteeheea hs 161 
Tables: Fy: (DIS OD): — xaes sn cabeeredersiacanceng cus ccaedte teshagea cecetathncs bunicese abana agence aseneag oan es 162 
Table6e:FQGsb). wh date dnote te Gast n Side ntalece hs Gain a barta aes 163 
Pable:PJ- «st san Sad bbvacetanrian ol aviel eznanidndelsceascheaniinsncavcedatadeaiusadevemestecouzedannesuagecs 164 


Publisher’s Notes: Great care was taken to transcribe Dr. Shimizu’s handwritten tables; however, errors may 
still exist. If you find something in question, please feel free to contact the publisher, and we can re-check the 
original manuscript. Dr. Shimizu edited and read the text several times. He found several of my typographical 
errors. He also added a few notes where his thinking had changed over the years. 


The original photographs were not available, so copies of Xeroxed copies had to be used. 


A special acknowledgement goes to Barry Sturman, who spent many hours copy editing the text and proofing the 
multitude of formulas. His wife, Susan, was also very considerate to give up so much family time to this project. 


Page iv Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


Chapter 1 — Preface 


The beauty of fireworks arises from their 
creation of organized patterns, mainly of fire, 
that consist of colors, lines and their movement 
and change in space and time. Chrysanthemum 
shells, also called ‘warimono’ (which means 
‘hard-breaking shell’), produce uniquely beau- 
tiful effects by creating patterns resembling 
chrysanthemums on the background of the sky. 


Until now, the design of chrysanthemum 
shells has been based entirely on the manufac- 
turer’s experience and intuition. The purpose 
of this paper is to introduce a more scientific 
method, using experiments and calculations. 


It must be noted that the requirements for a 
beautiful fireworks effect are not always strictly 
physical, but also sensual. For example, physics 
tells us that the trajectory of a star can never be 
completely straight. All that the artist requires, 
however, is that it should look straight. 


The stars that produce light for fireworks 
effects are well developed, but the characteris- 
tics of the bursting charge are not as well un- 
derstood. One of the important aims of this 
paper is to explain how to select the most ap- 
propriate bursting charge. 


This paper is divided into seven chapters. 
Chapter 2 describes the basic facts about the 
construction of the chrysanthemum shell. Chap- 
ter 3 studies the problems of design. Chapter 4 
discusses the preparation for experiments and 
necessary pre-calculations. Chapter 5 describes 
the experiments of burning the bursting charge 
in atmospheric air as a supplement to Chap- 
ter 4. Chapter 6 derives the formulas for the 
initial and traveling velocities of stars in the air. 
Chapter 7 prepares ballistic tables for moving 
stars and explains the relationships of the fac- 
tors that affect the velocity of stars. Further, it 
describes the modification and safety aspects 
of the bursting charge. Chapter 8, using the for- 
mulas and ballistic tables, compares the ex- 
perimental and calculated ballistics of stars. 
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Chapter 2 — Construction of Chrysanthemum Shells 


Japanese firework shells are of two principal 
types, the so-called poka and the warimono 
(chrysanthemum). The poka breaks just hard 
enough to disperse its contents into the air. The 
shell is made only as strong as it needs to be to 
withstand the shock of being fired. In contrast, 
the chrysanthemum shell must be strong enough 
to allow the bursting charge to distribute the stars 
far from the point of explosion. The quantity of 
the bursting charge must be great enough to meet 
this requirement. There are, however, shells hav- 
ing characteristics between those of warimono 
and poka. 


A typical warimono is shown in Figure 1. It 
consists of four main parts: shell casing (A), stars 
(B), bursting charge (C), and time fuse (D). The 
shell casing holds the contents and enhances the 
effect of the bursting charge. The stars are ignited 
and fly radially, producing a chrysanthemum 
image. Initially, the bursting charge is ignited by 
the time fuse, then the expanding burning gas 
breaks the shell and ejects the stars. At the same 
time, it ignites the stars. The time fuse is ignited 
by the flame of the lift charge and provides the 
delay before the explosion of the shell after it is 
fired from the mortar. The construction pres- 
ently used is as follows: 


The shell casing consists of two layers: the 
outside (a) and the inside (b) shell. Often, only 
the inside shell is called “the shell” (tamakawa). 
The inside shell holds the contents and keeps 
the stars arranged to produce a chrysanthemum 
pattern when the finished shell explodes. Con- 
sequently, the inside shell does not need to be 
particularly strong. The best material for the 
inside shell is newspaper, which is layered to a 
certain thickness on a wooden ball, and then it 
is cut into two pieces (hemispheres). The thick- 
ness depends upon the size of the shell, but 
usually it is 2 to 3 mm. The role of the outside 
shell is to strengthen the effect of the bursting 
charge, and it must have great tensile strength. 
The thickness of the outside shell is obtained by 
pasting narrow strips of strong paper over and 
over again on the inside shell to obtain the de- 
sired strength. In the old times in Japan they 


Figure 1. Section of a chrysanthemum shell. 


used fairly thick Japanese paper; however, at 
present, Kraft paper is used. An example of the 
pasting operation is as follows: 


Firstly, to avoid separation of the inside shell, 
in which the bursting charge and stars are ar- 
ranged in a pattern, a long paper tape is pasted 
along the joint between the two hemispheres. 
We call the operation dobari. Then taking two 
poles —A and A— on each hemisphere as 
guides, the paper strips are pasted on the surface 
of the sphere like the stripes of a watermelon. 
(See Figure 2.) The paper strips are cut so that 
the direction of the fibers is at a right angle to 
the length of the strip so that it conforms well to 
the inside shell. To make the outside shell strong, 
the width of the paper strips, if possible, should 
be constant. In Figure 2, (A) is the side view and 
(B) the top view. The numbers on each strip show 
the order of pasting. Pastings 1, 2 and 3 are pasted 
at one time. Next, new poles are taken at posi- 
tions that are perpendicular to the first pole axes. 
After two or three pastings, the shell is placed 
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(A) 
Figure 2. Outside view of pasting on the shell. 


between two thick wooden plates, and rolled to 
remove any projections and make the surface 
smooth. The operation is called gorokake. After 
this operation, the shell is dried in the sun. When 
the shell is well dried, the same processes are 
repeated until the desired thickness of the out- 
side shell is obtained. The last pasting is finished 
very carefully for good appearance. The num- 
ber of layers of paper for each sun, using Japa- 
nese paper, is 8 or 10, based on the tradition of 
Japanese fireworkers. (Sun is a Japanese unit, | 
sun = 3.03 cm.) For example for a 5 sun shell 
(15 cm shell): ‘sun 8’ means 5x8 = 40 layers. 
This method (i.e., to increase the strength of the 
shell in proportion to the diameter of the shell) 
is theoretically correct, as described later. 


(B) 


Fireworkers put their greatest effort into 
making the stars, because they are the main com- 
ponent of the beauty of fireworks. The stars that 
are used for fireworks are divided into three 
types: pumped cylinder stars, cut stars, and round 
stars. (See Figure 3.) 


The pumped star is shaped using a cylinder 
and a piston. For this manufacturing process the 
star mixture is dampened with a small amount 
of water, using glutinous rice starch, etc. as a 
binder. Cut stars are manufactured by making a 
cake of composition on a wooden board that 
has edges of the desired height. A dampened 
mixture is placed between the edges and the 
composition is flattened by tapping with a 
wooden hammer. The cake is cut into cubes and 


Figure 3. Fireworks stars: (A) pumped cylinder star, (B) round star, and (C) cut star. 
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well dried. Cut stars are often used for warimono 
because the color-changing stars are manufac- 
tured by covering 3—5 mm cubes with some 
colored mixture to form large, spherical balls. 


The operation is as follows: the cut stars are 
placed in a container with a spherical bottom. 
Some dry mixture is added to water to make a 
gruel (it is called toro). First the toro is sprin- 
kled over the cut stars and then it is uniformly 
spread over the stars, and a small quantity of 
the dry mixture is spread on the stars, and the 
container is shaken. The stars are somewhat 
round shaped. This operation is repeated two or 
three times. The stars should be dried quickly in 
the hot sun to prevent moisture from soaking 
into the already dry part of the star. When the 
stars are dry, they are cooled in the shade. 
These operations are repeated until the stars 
reach the desired diameter. 


A sample star that was made using such a 
process is shown in Figure 4. This star burns 
with a color change from red to blue and then 
flash. In Figure 4, (A) is the core (cut star) that 
is covered with an easily ignitable mixture to 
become (B). If (A) is not difficult to ignite, the 
same mixture (A) is used. This operation is 
called tomokokake, which means to cover the 
stars with the same mixture. 


The stars are next covered with a mixture of 
dark prime, which is the outer layer in (C). This 
layer is called “changing relay”. It used to pro- 
duce a good visual effect that makes the color 
change clearer. Then (C) is covered with a blue 
mixture of the proper thickness as shown in (D). 
Then (D) is covered with changing relay as 
shown in (E). Then (E) is covered with a red mix- 
ture as shown in (F). Finally, to ensure a perfect 
ignition, three steps are taken. The first is to 
apply a layer of relatively slow burning black 
powder type mixture on (F), the second is to add 
a layer of medium burn rate, the last is to add 
layer of fast burn rate to complete the star (G). 


The bursting charge is shaped into grains to 
be easily ignited. The bursting powder is mixed 
with 3 to 4% glutinous rice starch and a little 
water is added. This gruel mixture is poured 
over some grains of 4-5 mm, for example cot- 
tonseeds, to obtain round grains with a kernel in 
the middle. If cottonseeds are used, the grains 
must be dead and well dried for at least one 


o oO @ 
0 © €) 


(A) (B) (©) (G) 


Figure 4. The order of growth of the round 
Star. 


year. (See Figure 5.) For small shells, rice chaff 
is often used for the kernels. The weight ratio of 
cottonseeds to powder is 1:1 to 1:5. 


Figure 5. Bursting charge grains with cotton- 
seeds as kernels. 


The fuse for chrysanthemum shells is hand- 
made. The handmade fuse has a greater time 
delay and higher cost than commercial fuse. 
However, it has very good ignition properties, 
ignites other materials well, and causes the few- 
est ‘black’ shells (‘blind’ or ‘dud’ shells). 
Therefore, fireworkers willingly use the hand- 
made fuse. The core of this fuse generally con- 
sists of Japanese paper tape to which black 
powder is applied and which is firmly wound 
on a hard stick. This core is shown as e in Fig- 
ure 1. When the cover is removed, the core 
holds its shape. The flame from this bare part of 
the fuse becomes long and large, which can 
ignite even the least ignitable mixture. Also in 
Figure 1, the part marked f is split down the 
middle to ensure ignition from the flame of the 
lift charge. With this fuse, there is no dropout 
of powder from the fuse. With commercial fuse, 
such a procedure is impossible. The manufac- 
turer can easily regulate the burn time of hand- 
made fuse by changing the composition of the 
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Figure 6. An example of assembling a chrysanthemum shell. 


fuse mixture as required. This is not easy with 
commercial fuse. 


The fuse of the chrysanthemum must be 
longer than a certain defined length, adapted to 
the size of the shell. Commercial fuse, there- 
fore, is not often acceptable. 


Handmade fuse can be made as follows: Us- 
ing a brush, both sides of a sheet of Japanese 
paper are pasted with dampened black powder 
gruel and then dried in the sun. It is called ya- 
kushi (powdered paper). This is cut into narrow 
paper strips, about 80 mm wide and 250-300 
mm long, both sides are spread uniformly with 
powdered black powder. Then a thin bamboo 
stick is attached to the edge of one side and the 
paper is rolled around the stick. The powder 
becomes the core. Then a sheet of kozo paper of 
the same length as the core and of a proper 
width is attached to the core and rolled on the 
core. Then on a worktable it is slightly tight- 
ened by hand. Then the bamboo stick is re- 
moved. The kozo paper becomes the core cover. 
The paper stick with the core powder is then 
placed on a worktable and, using a thick wooden 
plate, the paper stick is pressed and rolled in 


one direction repeatedly. This is continued until 
the stick is rolled firmly. The worker should per- 
form about 100 of these rolling operations on a 
fuse that is 250 mm long. The first rolling op- 
eration is finished in this way. The second roll- 
ing process should be continued as it was for 
the first rolling. The fuse manufactured in this 
way burns 1.0—-1.2 seconds per | cm length. 


An example of assembling an ordinary chry- 
santhemum shell is shown in Figure 6. First the 
inside shell is punched at the center, and a piece 
of time fuse is inserted, and it is attached with a 
piece of hemp twine (d). Then a short paper 
tape is pasted over the outside end of the fuse to 
protect the end of the fuse. In the corner of the 
paper, the name of the shell or its contents is 
recorded. On the interior end of the time fuse, a 
paper bag (£) is attached that contains the 
bursting charge. The outer bag (A) is pasted on 
the inside of the shell casing, and it is turned 
inside out. This is used to hold the stars of the 
upper hemisphere. 


For the bursting charge bag (E) and the out- 
side bag (k) strong Japanese paper is used. 
These are placed on a stand (g) as in (A); (B) 
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shows the stars arranged along the inside of the 
lower hemisphere. The appropriate weight of 
bursting charge is then placed in the bursting 
charge bag and it is tapped on the outside with a 
piece of wood until there is no clearance be- 
tween the bag and the shell. The top of the bag 
is closed with hemp twine at (i), forming the 
upper part of the bursting charge bag that is 
spherical in shape. Then the bag (A) is turned 
up. The upper stars are arranged with the help 


of the bag (A) as shown in (D). The upper inside 
hemisphere is put on as shown in (E), and the 
upper hemisphere is tapped with a piece of 
wood to settle the position of the stars until 
both hemispheres make contact with each other 
as shown in (F). The marks on each hemisphere 
coincide with each other. The assembled shell 
is then pasted with paper tape to protect it from 
falling apart, and it is sent to the next operation 
of pasting the outside with paper strips. 
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Chapter 3 — Problems in Designing Chrysanthemum Shells 


In manufacturing chrysanthemum shells, 
most of the effort of fireworkers has been di- 
rected towards obtaining a certain number of 
‘petals’ and at the same time obtaining a large 
flower radius. In this paper the author evaluates 
the fundamental conditions required to meet 
these objectives. 


The problems are two-fold. The first prob- 
lem is the initial velocity attained by the stars 
from the action of the burst charge, that is, the 
problem deals with the explosion (or bursting) 
of the shell. The next problem is the motion of 
the star, which, once initially accelerated, flies 
in air with a trajectory of some sort and travels 
a certain distance. In other words, the problem 
is one of ballistics. 


3.1. Factors that Have an Influence 
on the Initial Velocity of Stars 


It is difficult to judge the influence of the 
various factors without experiments. However, 
from qualitative considerations they are thought 
to be as follows. The burst charge begins to burn 
upon its ignition. When the pressure reaches the 
breaking pressure, the shell explodes. The stars 
begin to move as a result of the action of the 
expanding gas from the explosion. This accel- 
eration ceases when the pressure behind the star 
becomes equal to the pressure that comes from 
the air resistance to the motion of star. It is 
known from interior ballistics’! that the initial 
velocity of the star is affected by the following 
factors: (1) the maximum pressure of breaking 
the shell, (2) the loading density of the burst 
charge, (3) the sectional density of the star (the 
value of the mass of star divided by the star’s 
largest cross-sectional area), (4) the force of the 
explosion of the burst charge and (5) the burn 
rate of burst charge at normal air pressure (1.e., 
vivacity). 


Generally, the combustion reaction of the 
burst charge is most effective in producing a 
high pressure when the breaking strength of the 
shell is large. This imparts a high initial veloc- 


ity to the stars. It is necessary to increase the 
strength of the shell for the stars to travel a 
greater distance. 


3.2. Calculation of the Break 
Strength of a Shell 


The breaking resistance of a hollow spherical 
body (such as a firework shell) is introduced as 
follows. 


Symbols are defined as: 


r; = inner radius P; = internal pressure 


r,= outer radius P,,= external pressure 


r = intermediate radius between r; and r, 


On the shell, the internal pressure P; and the 
external pressure P, act on each other at the 
same time. The internal pressure P; is a force on 
a unit area of the inside surface acting along the 
radius and is uniformly distributed on the sur- 
face. P, is the external force that acts on a unit 
area of surface along the radius of the sphere 
and is uniformly distributed on the outside sur- 
face of the shell (Figure 7). 


Figure 7. The pressures that act on the paper 
shell. 
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The following two assumptions are made 
about what happens when the sphere of the shell 
is deformed. 


(1) When the sphere is deformed by the ac- 
tion of a force from the outside or inside, the 
molecules on either surface of the sphere re- 
main on that same surface. 


(2) Consequently, after the deformation, 
the molecules originally on a surface of radius r 
exist on the deformed surface of radius r + Ar 


From above assumptions, the smallest vol- 
ume dv in the wall of the paper shell is assumed 
to be a small spherical hexahedron as indicated 
by ABCD-EFGH in Figure 8. 


Center of 
i the Shell 


Figure 8. A small spherical hexahedron 
representing the thickness of the paper shell. 


In this case 
AB = rdg = EF 
AE = rd@ = BF 

Therefore, the volume of the hexahedron is 
dv = r'drd@do 


When no pressure acts on the surface, the 
length AO and AD are r and dr. When pressure 
acts on the surface, these values change to 
r+ Ar and dr + d(Ar) [i.e., the both lengths in- 


crease: AO increases by Ar, and AD increases 
by d(Ar)]. Therefore, the deforming ratio of the 
molecules along the radius is d(Ar)/dg. When 
no pressure acts, the length of AE is rd@, and 
when the pressure acts it becomes (r +Ar)d0 
and it deforms to Ard@. Therefore, the deforma- 
tion in the direction of the tangent is Ard@/rd0 = 
Arir. 


Three surfaces at the point A are denoted as: 


ABFE = r’dgd0 = f, 
ADHE = rdrd0 = f, (1) 
ABCD = rdrdg = f, 


The surfaces of the small volume dv meet at 
right angles to each other before the deforma- 
tion. After the deformation the surfaces are at 
right angles to each other as before. Therefore 
the volume is not affected by shearing force, 
but only by internal force at a right angle. The 
internal forces received by the surfaces f,, f, and 
f- are denoted as 6,, 6,, and ©-. 


The small volume dv changes by the action 
of 0,, 6,, and 6, so that the unit lengths are de- 
formed, becoming d(Ar)/dr, dr/r, dr/r. Then the 
internal forces are calculated as: 


(1+m)—— eee ia + 
= fc =e dr r 
*  m(3-m) Ge 
r 
ae: m0 m=" + 
*  m(3-m) ees 
r 
That is, 
eee d(Ar) 
OSS Ga {c+m se +2(1—m) = 
(2) 
SS goke = 2 d(ar) 
9, =0.=— I ~) dr 2} ay 


Here E is the elasticity, and m is the value 
for the volume deformation derived from the 
deformation coefficient of the length along the 
direction of the outside pressure. Following the 


Page 8 Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


above equations, the interior force has no rela- 
tion to @ or 8. 


The values of the deformation coefficient, 
d(Ar)dr, Ar/r, of the volume dv along the three 
axes are then determined. For this purpose, first 
an equilibrium formula of the internal forces on 
the micro cube is constructed. The interior force 
on the three surfaces are 6, ft, 06, fi, O- f. The 
area of DCGH, which faces ABFE, is 


f. =(r+dvy dgde 
and the force acting on the right angle is 
o,, =o, +(dAr/dr)dr 


Therefore the total interior force acting on: 


oe \«, (S22) arf 
dr 
= \c, (S22) arp dr)? 


The surface BCGF is equal to the area of ADHE, 
f, = fy =rdrde@ 


(3) 


and the interior right angle force 6, is equal to 
the right angle force 6, on the surface of ADHE. 
The value of 0, from equation 2 is independent 
of 6 and a function only of 7 and the surface of 
radius r, the values of 6, should be equal in 
every point. Therefore, it becomes 


o,f, =O, 
In the same is true for the surface EFGH: 
O.f. = O.f, 


Therefore on the surface of the micro cube 
there arise six forces at interior right angle with 
no shearing force. These forces should balance 
each other as seen in Figure 9. 


Figure 9 shows the surface ABCD, but the 
other surfaces are in the same state. For the 
equilibrium of these forces the following rules 


apply: 


wa 
“ 


O 


Figure 9. Equilibrium of internal forces in the 
thickness of the paper shell. 


(1) When the forces are projected onto the 
three right angle axes, the sum of the compo- 
nents of force is zero. 


(2) The sum of moments along the three 
right angle axes is zero. 


The equilibrium equations of internal forces 
are as follows: 


for the direction of X axis: 
LX = On f. = Cfy a. Of, sin(4d¢@) 
- Of, sin(4dg) — o,f, sin(4dé) 


(4) 
o..f,,sin(4dé) 
=0 
for the direction of Y axis: 
LY =0,,f,,cos(4dg)—o, f, cos(4Zdg) (5) 


=0 


the direction of the Z-axis should be the same 
as for the Y-axis. 


From equation 5 we have o,,f,, =o, f, . This 


is the relationship described earlier. In equation 4 
relationships (1) and (3) are substituted, and we 
have 
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2(0, 25g. s 
r 


In this equation the multiples of dr can be omit- 
ted to give: 


Ac, -o,) +780 =i (6) 
dr 


Equation 2 is differentiated with respect to r 
and multiplied on both sides by r: 


2 
(+m) 4 


jf SO oa = (7) 
dr m(3-m) dr 


20 - Pool 
y 


When equations 2, 2' and 7 are incorporated 
into equation 6, we have 


2 
d (Ar) 15. 1 d(Ar) Ar 


dr? r dr r = (8) 
Now we put a special solution for equation 8 as 
Ar=r" (9) 
and the value of 7 is obtained: 
SOO = nr p FO) nny (10) 


These are introduced into (8) and we have 
n(n—l)r"? +2nr”? —2r"? =0 
The value of 7 is not zero, and it must be 
n(n—-1l)+2n-2 = 0 
Therefore 
(n—-1)(n+2)=0 


Solving this, we have 


ny = 1 (11) 
or 
ny =—2 (12) 


When these are introduced into equation 9, we 
have two linearly independent special solutions: 
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Arn=r™ , An=r” 


Therefore, two constants A and B are intro- 
duced and we have a general solution of equa- 
tion 8: 


Ar = Ar" +Br'? (13) 


Next the formula to obtain the equivalent 
stress is established. First, the constants A and B 
are obtained as follows: 


From equation 2 


o. =ef(l+m a +2(1—m) ~ (14) 
r r 


where 


E 


ae) (15) 


Equation 13 is differentiated with respect to r 


d(Ar) _ 


=n,Ar™ 7! + Br”! (16) 
dr 
Equation 13 is divided by r 
AP = Ap"! + Byte (17) 
} 


Equations 16 and 17 are introduced into equa- 
tion 14 


O,= e{(I +m)nr"'A+(1+ m)n,r™' BI 
+2(1—m)r"'A+2(1—m)r™"B (18) 
=d (Aar™" +BBr?" ) 
In equation 18 


fe eee 


(19) 
B=(1+m)n, + 2(1—m) 

The value of the normal internal force is equal 
to the internal force on the internal surface and 
to the outside pressure P, on the outside of shell 
with. Both of the pressures on the inside and 
outside of the shell press the thickness of the 
shell along the direction of the radius, and when 
r =F; O= —P; and when r = r,, and 0, = —P,. 
Therefore from equation 18 


-P= e( ar" +BBr?" ) (20) 
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-P, =e(Aar,"' + BBr,"") (21) 


Equation 20 is multiplied by 77”, equation 21 
is multiplied by 7,7,” , and then they are sub- 
tracted from each other: 


Ll 


Pr, Zn za Pre = eBB ( ge r = ae re ) 
From the above: 


my nm 
pow fatal Fin 
eB n. a — yp? rn 


a L 


(22) 


With the same calculation, multiplying equa- 
tion 20 by rz, and equation 21 by 7,7," and 
subtracting from each other we have: 


Pig Pre ea Aer? 17") 
From the above: 
Nn nN 
A= 1 bat og é a ae (23) 


When equations 22 and 23 are introduced into 
equation 17 we have: 


(P rr —Prrs ) (24) 


Next, when the shell is broken, it is suffi- 
cient to consider the deformation resistance 
along the direction of the tangent. This force is 
denoted as S,,. 


Table 1. Calculated Values of S, /P.. 
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Ar 
S,=E-— 
r 

Ny 

(Prn® — Pre)" 

_  m(3-m) B 

my Maly my -1 
Beg tg 


(25) 


The force comes only from the inner pressure. 
Therefore P, = 0. 


pesca [Sao] (26) 


Ji NM __ 4,224, 
Be Ny Ag 


An example of the calculation is as follows: 
m= 1/3 


(a) 7; = 68 mm, r, = 76 mm (corresponds to 
chrysanthemum No. 5), and 

(b) r; = 53 mm, 7, = 61 mm (corresponds to 
chrysanthemum No. 4). The values of S, /P 


are calculated, and the results are shown in Ta- 
ble 1 and Figure 10. 


50 60 70 80 
ymm 
Figure 10. Calculated values of S, /P.. 
(See Table 1). 


From the results of the calculations in Ta- 
ble 1 and Figure 10, we see that the resistance 
to the deformation is the greatest at the interior 
surface. Therefore, when the value S a is greater 


than the internal resistance force, the shell is 
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Table 2. The Critical Internal Breaking Pressure P, (kg/cm’) for Internal Diameter y; and 


Thickness of Shell 6. 


broken. Therefore, we can obtain the critical 
internal breaking pressure when we set r = 7;, 
S,, = Ji. The critical limit of the internal break- 


ing pressure is given by equation 27. 


oe eg oF Te (27) 
= 
i Slee Os ie OL 
a Bp 


Example: J, = 210 kg/cm’ (internal force of 
resistance of Japanese paper), the internal ra- 
dius of the shell 7; = 25 — 150 mm, the thickness 
of the shell = 2 — 10 mm: the calculation results 


60 


Figure 11. Graphs for Table 2. 
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are shown in Table 2 and Figure 11. Then, for 
various internal radii the thickness and the dif- 
ference of the thickness for a definite breaking 
force are shown in Table 3. 


According to Table 3 the values of Aéd are 
constant for a definite interior pressure, even 
when the values of 7; change. Consequently, 
when the internal breaking pressure is in- 
creased, the strength of the shell should be in- 
creased in proportion to the radius of the shell. 
This coincides with the experience of firework- 
ers described in Chapter 2. 


100 120 140 160 
P, (kg/cm?) 
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Table 3. Internal Diameter y; and Thickness of Shell 5 and Aé for Critical Pressure y¥,. 


| 1 20 50 


(mm) 5 | Aa | 8 


Py (kg/cm?) 


25 


It should be noted that equation 27 applies 
to a theoretically uniform material that forms a 
hollow sphere. However, shells are not uniform, 
because they are made of paper. The calcula- 
tions just described must be expected to have 
some unknown problems when they are applied 
to the design of a shell. The value for “m” 
(Poisson’s number) should be determined ex- 
perimentally.”!_ None the less, this method 
should be useful as a rough guide to designing a 
chrysanthemum shell. 


3.3. The Effect of the Type of 
Paper on the Strength of the 
Shell 


Experience teaches us nothing except the 
principle of per sun, which means the strength 
of the shell is proportional to the folds of paper 
pasted on the shell without consideration to the 
diameter of the shell. That is to say, for exam- 
ple, eight folds per sun of Kraft paper on the 5 
sun (6 inch) shell gives the same breaking 
strength as eight folds per sun on the 10-sun 
(12-inch) shell. However, it does not teach us 
about the quality of the paper. 


In Japan we use Japanese paper (Kozo Pa- 
per) or Kraft paper. Japanese paper is soft and 
before use, two or three sheets are pasted to- 
gether into one sheet. It is easy to work with, 
but its high price prevents using it. Kraft paper 
is less convenient to use than the Japanese pa- 
per, but it is less expensive, so at present many 
fireworkers use Kraft. In this paper both types 
of paper are studied. 
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3.4. Mass per Unit Cross Sectional 
Area of a Spherical Star 
(Sectional Density) 


The equation for the start of a star’s motion 
is as follows: 


d’x 


m=*=0(P,-P,) (3) 


where m is the mass of a star, o is the area of 
the greatest cross section of the star, d’x/d? is 
the acceleration, P;is the resistant force per unit 
area (pressure) that acts on the front of the star, 
and P, is the pushing pressure that acts on the 
rear of the star. 


Integrating equation 28 we have: 
av=" "(BR -P) dt, Av=v—v,=0 (29) 
m-? 9 


This shows that the change in the velocity of 
the star is inversely proportional to the value of 
the sectional density m/o, which depends on 
the true density and the size of the star. When 
the true density of the star is constant, the value 
of m/o increases or decreases in proportion to 
the diameter of the star. The true density is gen- 
erally 1.3-1.6 grams per cubic centimeter, and 
it is difficult to obtain larger or smaller values. 
The smaller the sectional density, the larger the 
initial velocity of the star. However, after the 
driving force of the burst charge diminishes, the 
star is greatly decelerated by the air resistance. 
This must be remembered when designing the 
chrysanthemum. 


3.5. The Compositions of the Burst 
Charge 


The compositions that have long been used 
in Japan were mixtures of potassium chlorate 
and hemp charcoal, which were called the hai- 
san explosives. (The meaning of the word hai- 
san or H3 comes from the mixing ratio of hemp 
charcoal 3 parts to potassium chlorate 10 parts). 
For safety reasons some /aisan are mixed with 
potassium nitrate. It must be remembered not to 
mix sulfur or red phosphorus with the haisan. 
In ordinary firework factories, they use Black 
Powder that contains sulfur. In spite of the dan- 


Page 14 


ger of mixing sulfur with chlorate, the reasons 
for using haisan are: 


1) it is very quick burning, 
2) it provides a large explosive force, 


3) it requires only a simple manufacturing proc- 
ess, when compared to black powder, and 


4) it is the cheapest of all. 


In this study three compositions are studied: 
1) Black Powder, 
2) potassium perchlorate composition, and 


3) potassium chlorate composition (haisan). 


3.6. The Force of Explosives f 


Following the theory of internal ballistics, 
the energy output E of a propellant is denoted 
as 


pe (30) 
y-1 
where w is the mass of the explosive charge, fis 
the energy output per unit mass of explosive 
and 7 is the average adiabatic expansion coeffi- 
cient of the gas, typically 1.25. 


The value of fin joules/kilogram is given by 
nRT where n is the number of moles of gas pro- 
duced by | kilogram of explosive, R is the gas 
constant and T is the explosion temperature in 
kelvin. If Vp is the volume of gas produced by 
the explosion of 1 kg of explosive, reduced to 
273.15 K and | atmosphere pressure, and if the 
energy is expressed as kilograms weight x 
decimeters, the value of fis given by 0.3782 x 
V, x T, units decimeters. This is how f will be 
calculated in Chapter 4. 


When burnt gas of the burst charge projects 
a star of mass m with the initial velocity V we 
have: 


¥ymV* = ey oR 
y-1 
Here, R denotes the effective fraction of the 
charge that acted on the unit cross sectional 
area 6 of the star, and it is called the efficiency. 
From the above equation: 
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(31) 


1 H “hk 
oO 


The value of fw is independent of R, and 


the initial velocity is proportional to the square 
root of fw. This should be confirmed by ex- 


periments. 


3.7. The Vivacity A of Burst Charge 


The vivacity of the burst charge has an in- 
fluence on the initial velocity of a star with the 
effect of the breaking strength of the shell. Us- 
able quantitative data concerning this problem 
have not yet been obtained. The core substance 
of the burst charge should be also selected to 
give good values for the vivacity. This is the 
important goal of this chapter. 


3.8. Loading Density 


The construction of the burst charge of the 
chrysanthemum is almost uniform among fire- 
workers, and the values are 0.28-0.32 kg/dm’. 
The values are obtained from the mass of the 
burst charge divided by the space occupied by 
the stars and the volume of supporting material 
(core). The relation of the breaking strength and 
the loading density is denoted by the equation 
of Noble and Abel:!"! 


Rafe ho (32) 
ae 
A 


where A is loading density, Z) is the burning 
ratio of the burst charge at the breaking of the 
shell, namely the ratio of burnt burst charge to 
initial mass of the burst charge, n the covolume. 
The breaking pressure Py is thought to be con- 
stant, and Zy decreases with the increase of A. 


3.9. Factors Concerning the Velocity 
of the Stars 


Air decelerates moving stars. According to 
fluid mechanics, air resistance acting on a star 
is given by: 
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6 6 
D=K' pvo[ -k'po( 2 (33) 
ply Re 


where p is the density of the air, p is the coeffi- 
cient of viscosity of air, v is the velocity of the 
star, ois the cross sectional area of the star, @ 
is the diameter of the star, and Re is the Rey- 
nolds number: 

Re= ply at 0 ae Force of inertia 


eg 


Force of viscosity 
(34) 

The fluid is air. Therefore, the force of inertia is 

much greater than the force of viscosity. The 


value of the Reynolds number is therefore also 
very large. From equation 33 


Be 


From this the value of 6 is obtained: 


eas (35) 
log Re 


As the value of Re approaches infinity, the 
value of 6 approaches zero. 


The terms 


Gis) ~#(a) 


can be replaced by another constant, k. 


The motion equation for a star flying in a hori- 
zontal direction is therefore 


(36) 


For a star that flies in a vertical direction the 
motion equation is: 
d 
m— =-kuo — mg (37) 
dt 
where g is the acceleration due to gravity. The 
above two equations are summarized as fol- 
lows: 
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dv ky 


dt pit) (38) 
ee 
dt p(t) = 


where p’ is the sectional density of the star, i.e 
the mass of the star divided by the cross sec- 
tional area of the sphere at that time. 


3.10. The Relation of Physical 
Conditions and the Visual 
Beauty 


As described earlier, the uniformity and the 
straightness of the motion of star are most im- 
portant for the beauty of the chrysanthemum. 
The trajectory of stars should be as straight as 
possible and the trace of the moving stars 
should be uniform with respect to time and 
space. 


The chrysanthemum must develop the 
beauty of the star light as a quantitative effect. 
Therefore the number of stars contained in a 
shell should be as large as possible. If the di- 
ameter of the shell is limited, the stars should 
be as small as possible. However, in this case 
the density of the stars becomes smaller as the 
diameter of the chrysanthemum becomes 
smaller. These factors should be examined later. 


3.11. Calculating the Number of 
Stars in a Round Shell 


The number of round stars that should be 
contained in a spherical shell is calculated as 
follows. The diameter of the shell is denoted as 
D, and the diameter of the stars as d. Now six 
stars are arranged in a flat layer as shown in 
Figure 12. Actually, there are seven stars, the 
centers of which are o, a, b, c, d, e, and f. 
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Figure 12. Basic calculation of the number of 
stars arranged in a shell. 


The area of the hexagon abcdef is: 


V3 


S=—d’x6 
4 


This area contains the center star “o” and 1/3 of 
a hexagon ABCDEF. Therefore this area is the 
space that can contain 3 stars (1 + (1/3) x 6 = 3) 
and the area per star is: 

3 

G=4S= 3p 

2 
However, the centers of the stars are not on a 
flat surface, but on a sphere, so the above for- 
mulas only show the rough calculation. 


The distance from the center of shell to the 
contact point of each star is “2(D-d) cos 9, 
where 6 is the angle of a star from the center of 
the shell. The area of the sphere that contains 
the contact points of the stars is: 


6’ =m(D-d) cos’ @ 
= m(D-—d) (1-sin’ 0) 


=2(D-d) t a 


Therefore, the number of stars n is denoted 


as 
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Table 4. The Number of Main Stars along the Inside of a Shell. 


Internal Diam. Diameter of Number of 
Shell Size of Shell Star Stars 
(inch) (mm) (mm) (calculated) 


Number of 
Stars 
(empirical) 


ope Oe n=4352(2-2) (39) 
5 V3d?/2 
ap 2 DD This is a formula to calculate the number of 
= Sen (2 spherical stars to be contained in a shell. A 
comparison between the calculated and empiri- 
where & is an adjustment factor based on ex- cal number of stars are shown in Table 4. 


perience. If k is set to 1.20, we have: 
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Chapter 4 — Method for Experiments on Shell Burst, 
Preparation and Preliminary Calculations 


4.1. Experimental Method 


Three types of bursting charge were pre- 
pared: 


a) potassium perchlorate 
b) black powder 
c) potassium chlorate 


The stars were made with potassium perchlorate 
as the oxidizer. 


Sample 5- and 6-inch chrysanthemum shells 
were prepared with these materials. The shells 
were reinforced by pasting them with Japanese 
or Kraft paper. The samples were ranked based 
on the number of paper layers. 


Samples prepared in this way were exploded 
at a height of 1 m above the ground at night and 
the movement of the stars was photographed 
with a still camera (through a revolving shutter 
driven by a synchronous motor) and also with a 
16-mm movie camera. 


4.2. Aims of the Experiments 


(1) To obtain the relationships between the 
size of the star, density of the star, type of burst 
charge and the break strength of the shell. 


(2) To find the relationships between the ve- 
locity of the stars and air resistance, including 
the effect of the diameter, density, and burning 
velocity of a star. 
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4.3. Types of Sample 
Chrysanthemum Shells, their 
Construction and Preparation 

4.3.1. Stars 

(1) Types of Stars and their Symbols 


Six types of stars were prepared. They were 
identified with the following symbols: 


Dia.(mm) | Amber | Blue | Green 
Rall 8 Ce | Ba | Ge 
| 12 C 12 Bio G 12 | 


(2) Construction of Stars 


Using 3-mm flash stars as the core, compo- 
sitions SKC, SKB or SKG were pasted in the 
usual way for round stars. The flash core was 
used to show clearly the end point of the burn- 
ing. 
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(3) Compositions of the Stars 


SKCg (Amber Star) 


Potassium nitrate 


Glutinous rice starch 

SKB (Blue star) 

Potassium perchlorate 64.0% 
Pine root pitch | 13.0% 
Glutinous rice starch 6.0% 


SKG (Green star) 


| Potassium perchlorate 


Glutinous rice starch 6.0% 
SBF (Flash) 

Barium nitrate 67.0% 
‘Aluminum (fine flake) || 27.0% 
[Glutinous rice starch | 6.0% 


SBKF 


(pasted on the flash core to ignite easily) 


Barium nitrate 


Glutinous rice starch 

SKCg (Amber composition called 

Eight Chrysanthemum) 

Potassium nitrate 48.5% 


Glutinous rice starch 


Glutinous rice starch 6.9% 
SKCs (Ignited with Black Powder) 
Potassium nitrate 75.0% 


1.5% 
(add.%) 


Eight Chrysanthemum means the ratio of the char- 


coal to the oxidizer is 8:10 in weig 
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ht ratio. 


(4) Diameter of Stars 


Twenty five sample stars were selected from 
each manufactured batch of stars. The maxi- 
mum and minimum diameters of each star were 
measured. In the following tables the maximum 
and minimum, or minimum and maximum, di- 
ameters are recorded for each sample star. 


Cs (Amber star, 8 mm dia.) 


The difference between the measured value 
and the average is denoted by v, and the number 
of the stars is denoted by n. The probability 
deviation r of the diameters of the stars is calcu- 
lated as follows: 


r =0.6745 
n— 
= 0.6745 0.048134 
50-1 
=0.15 mm 


Editor’s note: There is a 50% probability that an 
individual result will be in the range given by (aver- 
age t probability deviation). 
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Cy. (Amber Star, 12 mm dia.) By (Blue Star, 12 mm dia.) 


Ave. = 12.26 mm 12] 12.16 | 12.15 | Ave. = 12.11 mm 


= 018mm 14] 12.65 [12.34] r= 0.22mm 


Ave. = 8.73 mm : ‘ Ave. = 12.23 mm 


r=0.18 mm : : r=0.19 mm 
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F; (Flash Star, 5 mm dia.) 


Note: in this case only ten stars were sampled. 


(5) Mass, Density, Sectional Density, and Bulk 
Density (Mass of 1 Liter of Stars) 


100 sample stars from a manufactured batch 
were selected and weighed. The result was di- 
vided by 100 to obtain the average mass of one 
star. The density of the layers was calculated as 
shown below. 


Setting m, as the mass of a star, m2 as the 
mass of a flash core, 7; as the radius of a star, 
and r, as the radius of a flash core; then 6, the 
density of the colored zone and 6, the density of 
the flash zone were calculated using the follow- 
ing equations: 


In the same way, the initial sectional densities 
(i.e., mass of the star divided by the maximum 
cross-sectional area) p; and p’, were calculated: 


bo : 
Pi =a P2 = 
rr, nr; 


guce 


Introducing the data to these equations, we have 
the results shown in Table 5. The bulk density 
5! indicates the total mass of the stars that will 
fill a 1-liter container. 


Table 5. — Mass, m,, Density, — , Sectional Density, p;, and Bulk — gues of the Stars. 


(6) The Burn Time for Stars in the Air when 
the Stars Are Not Moving 


The stars were placed on the ground and ig- 
nited. The burn time was measured with a stop- 
watch. 
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Ave. = 1.59s 


Note: The stars marked with * stopped burning be- 
fore the flash, and they were omitted from the calcu- 
lation. 
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Ave. =2.17S 


Ave. =1.95s | 


Note: The stars marked with * stopped burning be- 
fore the flash, and they were omitted from the calcu- Gr 


lation. 
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OPWOPOINIODO:IO:iBiWwiNia 


Ave. = 0.60 s 
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4.3.2. Burst Charge 
(1) Types and Symbols 


Three types were prepared, identified with the 
following symbols: 


Mixing Ratio 
Symbol |Component (%) 


Potassium chlorate 
H 


Glutinous rice starch 2.0 (Add’l.) 


Potassium perchlorate | 70.0 


Hemp charcoal 18.0 

Sulfur 12.0 
Glutinous rice starch 2.0 (Add’l.) 
Potassium nitrate 

Hemp charcoal 

Su 

Glutinous rice starch 


(2) Construction of the Burst Charge 


The H and P burst charges were made as 
follows: The component materials were mixed 
by sieving, some water was added, and the re- 
sulting mixture was kneaded. The burst charge 
was then coated on cottonseeds and dried in the 
sun. The S burst charge was made as follows: 
The mixture was placed in a ball mill contain- 
ing wooden balls and milled for 18 hours. Then 
the mixture was coated on cottonseeds and dried 
in the sun. The construction data are as follows. 
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pH | P| s 


Weight of the charge 
(in kg) per 1 kg 1.313 | 1.625 | 1.288 
cottonseeds 


powdered charge 

completed charge 
The thickness of the burst charge that was 
coated on the cottonseeds was measured as fol- 
lows: twenty completed grains were selected 
and they were measured to determine the short- 
est, middle and the longest diameters, d), d>, d3, 
of each ellipsoid, then the grains were washed 
to remove their coatings and the lengths of the 
cottonseeds, d,', dy’, d;', were again measured. 


The thickness was calculated by the following 
formulae. 


Pa Sit 
2 

d, —d, 

i= a 2 

y= a>" 


The thickness (in mm) was obtained from (¢, + 
t) + t3)/3. The results of the measurement are as 
follows: 


Page 23 


Page 24 Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


From the above measurements, we see that 
although the thicknesses are very irregular, the 
average values are similar, and the values range 
from 0.47 to 0.49 mm. 


(3) Sensitiveness Test by Drop Hammer 


The results of drop hammer tests for the 
burst charges, H, P, and S are shown in Fig- 
ure 13. For comparison the data from tests of 
picric acid are also shown. The mass of the 
hammer was 2 kg and the thickness of the tin- 
foil cover of the sample was 0.05 mm. From 
the results we see that the S burst charge is the 
most insensitive and H and P have almost the 
same sensitiveness. The problem of fireworks 
safety is very complicated and these data only 
serve as an indication of the relative sensitive- 
ness. 


(4) Calculation of Characteristic Values of 
Burst Charges 


1. Calorimetric Data of Solids and Gases 
(a) The Internal Energy of Gases 


The internal energies at various temperatures 
were calculated by means of statistical me- 
chanics, using the normal vibration numbers 
from spectral analysis with the standard vibra- 
tion number at 15 °C [=E}j5°c (T)/J kcal/mol)]. 
They are shown in Table 6. 


However, for KCI the author did not have 
any reliable thermochemical data, and the cal- 
culation was conducted as follows: 


15-790 °C 
(E) = 0.126 (790-15) +5.52 
=15.29 kcal 
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Picric Acid 
90F (1953, 6.19, 


Percentage of Explosion 
oa 


= 
Oo Oo 


0 10 20 30 40 50 6 


H (1954, 
10.23, 
20 °C) 


0 70 80 90 100 110 120 130 


Height of Hammer (cm) 


Figure 13. The results of drop hammer test of burst charges. 


790-1420 °C 
0.0126 + 0.0060 


E)= 
(B) = {9 

= 5.86 kcal 
Greater than 1420 °C 


E=2R(T-1420-273)x2x—— 
2 1000 


} (1420-790) 


= 5.958(T -1693)x10° 
=5.96x10°T —10.09 kcal 


The sum of the above should be the formula 
that is applicable for the temperature range up 
to 1420 °C: 


Exc = 11.06 + 5.96 x 10° 7 


(b) The Internal Energy of Solids 


E is the internal energy of 1 gram atom 
and C, is the atomic heat, then 


Here x= and © is called the specific tem- 


perature. 


Table 6. The Internal Energy of Gases with Temperatures E15 °c (T)/J kcal/mole . 
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Table 7. The Values of C, and £/© as a Function of 1/x = 7/0. 


The values calculated from the above for- 
mula of C, and £/© as a function of 1/x = 
T/© are shown in Table 7. 


In the following paragraphs, brackets [ ] in- 
dicate the solid state and parentheses ( ) indi- 
cate the liquid state for a substance at normal 
temperature. 


In the calculations, four substances have to 
be included: carbon [C], potassium sulfate 
[K,SO,], potassium oxide [KO], and potas- 
sium sulfide [KS]. 


When the atomic heat at high temperature is 
known, the value at a high temperature is di- 
vided by the number of atoms to obtain the 
value of C,. Then from C, in Table 7 the value 
of 7/© is obtained. T is known and © is cal- 
culated. To obtain the internal energy E£, the 
T/© value for each temperatures is calculated 


for ©. Then E/© is obtained from the col- 
umn of T/©. The internal energy £ is obtained 


by multiplying the value of E/O by © and 


the number of atoms in one molecule. The 
value of © is shown as follows.!?! 


Measured temp. T (K) 
Molecular specific heat 


Average atomic heat 
T/O 
© (kK) 


For [K,O] and [K,S] the value of © could 
not be obtained, so it was roughly calculated 
using the following formula. Using C, as 3nR, 
we have: 
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Ex,o = 0) = Ex, 3) = = 3nR(T - 288)x—— Table 8. The Internal Energy of Solid Ma- 


1000 terials E15 °c (kcal/mole ). 
1 


= 3x3x1.986% (I — 288)x 


=—5.154+17.87x10°T 


The values calculated in this way are shown in 
Table 8. 


2. The Characteristic Values of Bursting Charge (H) 
(a) The Composition of 1 kg Mixture Consisting of: 


Mixing Ratio | Quantity in 1 kg] Chemical | Composition | Quantity in 1 kg 
Materials (%) of Mixture (g) | Constituents (%) Mixture (g) 


Potassium chlorate 
Hemp charcoal 


C6H1905 


Glutinous rice starch 
HzO 


(b) Number of Moles of Components in 1 kg of mixture = x CO, + y H. +z CO 4 
1 kg Mixture u (HO) + a [KCl] + O 
KCIO; = 735/122.58 = 6.00 mole Here Q denotes combustion heat (explosion 
C = 201/12.01 = 16.74 mole heat). On the other hand, introducing the water 
H,O = 28/ 18.02 = 1.55 mole gas reaction: 
C6HpOs = 17/162.16 = 0.10 mole CO +H,0<_> CO; +H 
(c) Number of Gram Atoms in I kg Mixture Sue K(T) 
xy 


Using the following calculated data by statisti- 
cal mechanics:!°! where 


K(T)= CO* H,0 
EF CO, P H, 
[eco foof 2005 a7se[410} | 5. rans 


(d) Explosive Reaction 


The reaction is assumed to be as follows. 
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Table 9. Values of K(7) for Various Temperatures. 


(e ) Calculation of Characteristic Values 
from the Balance of Materials: 


C=x+z 

O =2x+zt+u 

H =2y+2u 
KCl =a 


From the heat of explosion Q 


Q= X40, + 24co + UG (u4,0) + Bde — 49 t Aq 
=XEco, + Ey, + ZEco + UE 44,0) + QE Key 
Here, g denotes the heat of formation; Aq, the 
value of the correction of g from constant pres- 


sure to constant volume and equal to 0.57 
kcal/mol gas (at the temperature of explosion). 


Ico, = 97.80 kcal 
Ico = 29.70 keal 
Yui,0) = 68.31 keal 
A«cy = 105.6 kcal 
Ikcio,) = 99-9 keal 
g = 6.00 dxcio, 1-55 Fai,0) + 9-10 Fo,.061 

= 6.00 95.9 +1.55 x 68.31+ 0.10 x 226.8 

= 704.0 kcal 


From the above formula we have: 
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y= C+ H-O+3x'=x-0.66 


z= C-x=17.34-x 
u= O- C-x=2.71-x 
Therefore 


Ag =0.57x(x+y+z+ut+a@) 
=0.5x{C+5H+a)=0.57%25.39 


=14.47 
0, = Xco, + 24co + U4 ,0) + CF Kc — 7 + Aq 
=97.80x + 29.70(17.34 - x) + 
68.31(2.71—x)+6.00x105.6 — 


704.0 +14.47 
=0.21x + 644.19 


Approximate Value of Explosion 
Temperature T,, 


When we disregard the mol of CO, (x = 0), 
we have Q; = 644 kcal. 


From the relation of internal energy, disre- 
garding x and y, we have: 


QO, = ZE oo + UEW.0 + OE gc 


Then we have: 


T= 2000 K Q;= 10.89 x 17.34 + 2.71 x 
(14.81 + 10.54) + 6 x 22.78 = 395.41 
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T=3000K Q,= 17.61 x 17.34 + 2.71 


ZE co UE Lo OE Ke Q2 
(T) | cay | @_ | eal) 


(25.97 + 10.54) + 6 x 28.94 = 7) | @ 


00059 OF | TC Ge 798 12810) 170 64 | 24859) 
3500 | 67-00 121-16 | 327.34 | 29.96 | 191.92 | 636.38) 


576.32 


T=4000 K QO) = 24.44 x 17.344 2.71 
x 
(36.54 + 10.54) + 6 x 34.90 = 
760.78 
( £y,o was added with the latent heat 10.54 
kcal/mol for vaporization.) 


When we read the point QO; = 644 kcal on 
the T vs O graph, we have an approximate value: 


T, = 3360 K 
More Accurate Calculation of T,,. 


From Table9 we have an equilibrium con- 
stant K(T) = 7.69 for 7, = 3360 K. Therefore: 


Zu (17.34—x)(2.71-x) -~769 
xy x(x-0.66) 


From the formula we have, x = 1.76, there- 
fore: 


y=x-—0.66=1.10 mol 
z=17.34-—x=15.58 mol 
u=2.71—x =0.95 mol 


Therefore: 
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‘4000 | 78.21 |24.77| 380.78 | 34.71 | 209.40 | 727.87 


O, =0.21x + 644.19 = 0.21X1.76 + 644.19 
= 645 kcal 


From Tables 6 and 8 having the values E and 
calculating Q2, we have: 


Reading the value for 7, at O2 = 645 kcal on 
the Tvs Q graph gives the explosion temperature 
T,. Thus we have: 7,,= 3550 K 


Repeating the calculation in the same way, 
using this value of 7,, the result did not change. 


The specific volume is calculated as: 
Vi =(xtytztu+a)x22.41 
=(1.76+1.10+15.58+ 0.95) 22.41 
= 569 | 
The force of explosives is therefore 
Fy =0.3782V.T, 
= 0.3782 x 5693550 
= (.736x10° dm 
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3. The Characteristic Values of Burst Charge (P) 
(a) The Composition 


Quantity in 
Materials 1 kg Mixture (g) 


Potassium perchlorate 


Glutinous rice starch 


(b) Number of Moles in 1 kg of Mixture 


KCIO, = 681/138.55 = 4.91 
C = 144/12.01 = 11.99 
S2 = 118/64.12 = 1.84 
H20 = 23/18.02 = 1.27 
CgHi00s = 17/162.14 = 0.10 


(c) Number of Gram Atoms in 1 kg Mixture 


(d) Explosive Reaction 


The reaction of the explosion would proceed 
according to the following process: 


Reaction (1) 


KCIO4 + S2 = KCl + 2 SO, 
1.84 1.84 1.84 3.68 
The remaining KC1O, is 


4.91 — 1.84 = 3.07 mole 
Reaction (2) 


KCIO, + 2C = KCl + 2CO, 
3.07 6.14 3.07 6.14 


Chemical | Composition Quantity in 
Constituents 1 kg Mixture (g) 
KCIO4 
H20 
others 


The remaining C is 
12.59 — 6.14 = 6.45 mole 
Reaction (3) 


Reduction of SO, 
SO, + 2C = 2CO + %S, 
3.23 6.45 6.45 1.62 


The remaining SO, is 
3.68 — 3.23 = 0.45 mole 
Reaction (4) 


CO, + SO, = CO + SO; 
0.45 0.45 0.45 0.45 


The remaining CO, is 
6.14 — 0.45 = 5.69 mole 


From the above calculations, the substances 
produced should be 0.45 mole SO3, 1.62 mole 
S, and 4.91 mole (1.84 + 3.07) KCI and CO, 
CO. The latter would be in equilibrium with 
H,O and Hb, in the gas reaction. Therefore, the 
reaction equation should be given as 


1 kg mixture =x CO. + y H,. +zCO+ 
u (H20) + 0.45 SO; + 1.62 (S2) + 
4.91 (KCl) + Q (reaction heat). 


From the number of atoms 


C=x+z=12.59 
O=2x+z2+u43x0.45 = 21.41 
H=2y+2u=3.54 
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From these relations 


y=x—-5.70 
z=12.59-x 
u=747-x 
(e) Calculation of Characteristic Values 


From the heat of formation, 
Q, = X40, + 24co FY ,0) F 0.45450, 
+4.9]ldye —g + Ag 


The relation of the equilibrium constant and the 
temperature comes from Table 9. 


First the author obtained 7, = 3900 K and using 
the corresponding value of K(T) = 8.36, a pre- 
cise calculation was carried out. 


z-u _(12.59-x)(7.47-x) 


= 8.36 
xy x(x-5.70) 


K(T)= 


From this equation we 
have x = 5.91, and then we 
have 


y= 5.91 -4.80=1.11 
z= 12,595.91 =6.68 
u=6.57—5.91 = 0.66 


On the other hand, 
Ico, = 97.80 kcal 
so, = 91.9 kcal 

GIco = 29.7 kcal 
A«cy = 105.6 kcal 
Yu,o) = 68.31 keal 


q= 4.91 dixcio,| + 1.274 4,0) + 0.10416.41,,0.] 
= 4.91x113.54+1.27x 68.31+ 
0.10 x 226.8 
= 666.72 kcal 


0.454+1.62+4.91 
5.914+1.114+6.68+ 
=0.57x 
0.66 + 0.45 +1.62+4.91 
=12.16 kcal 


xtyt+z+ut 
Aq =0.57x 


Therefore the value of Q; is: 
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Q; =5.91 x 97.80 + 6.68 + 24.70 + 
0.66 x 68.31 + 0.45 x 91.9 + 
4.91 x 105.6 — 666.72 + 12.16 


= 726.78 kcal 


From the internal energy of the component 
materials, we have: 


Q, = XE co, (T) + yEy, (L) + Zoo (F) + 
UE, o(T) + 0.45E go, (T) + 
1.62E, (T)+4.91E xq (7) 


From Table 6 the value of Q; was calculated for 
temperatures from 3000 to 4000 K. The case of 
E\,,o Was considered from the liquid state, and 


the value of the latent heat of evaporation of 
water 10.54 kcal/mole (20 °C) was added: 


A graph of T vs Q) was drawn and on the 
curve of T vs Q> the point corresponding to Q) 
= 727 showed the explosive temperature: T, = 
3940 K. 


The specific volume is calculated as: 


Vo =(x+y+z+ut+0.4541.62+4.91)x 22.41 
5.91+1.114+ 6.68 + 
7 ee 0.45 +1.62+4.91 
= 478 liter 


} 22.41 


The force of explosive is calculated as: 


fr =0.3782V.T, = 0.3782 x 478 x 3940 
=0.712x10° dm 
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4. The Characteristic Values of Black Powder (S) 


(a) The Composition of I kg Mixture: 


Hemp charcoal 


Sulfur 
Glutinous rice 
starch 


(b) Number of Moles in 1 kg Mixture 


KNO; = 735/101.10 = 7.27 mole 
C = 120/12.01 = 9.99 mole 

S> = 98/64.12 = 1.53 mole 

H,O = 18/18.02 = 1.00 mole 
Cs5H,pOs; = 17/162.16 = 0.10 mole 


(c) Number of Gram Atoms in I kg Mixture 


— | 0.50} 0.60} 1.00 


(d) Explosive Reaction 


As a hypothesis, it is proposed that the steps 
in the explosive reaction are as follows: [The 
number under the chemical formula shows the 
moles of the substance. | 


Reaction (1) 


2 KNO;3 + % So = K,SO, a N> + O, 
6.12 1.53 3.06 3.06 3.06 


The remaining KNO; is 
7.27 — 6.12 = 1.15 mole 


Reaction (2) 


2 KNO; + C = K,CO;3 + No + 3/2 Op 
1.15 0.58 0.58 0.58 0.87 


The remaining C is 
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10.59 — 0.58 = 10.01 mole 
The remaining O) is 
3.06 + 0.87 = 3.93 mole 
Reaction (3) 
C+ O, = CO) 
3.93 3.93 3.93 
The remaining C is 
10.01 — 3.93 = 6.08 mole 
Reaction (4) 


K,SO4 + 7/4 C = 4 KyCO3 + 4 KoS + % So + '%4 CO2 


3.06 5.36 1.53 153 0.77 


The remaining C is 
6.08 — 5.36 = 0.72 mole 
The remaining CO, is 
3.93 + 3.83 = 7.76 mole 
Reaction (5) 
CO, + C = 2CO 
0.72 0.72 1.44 
Thus the carbon is totally consumed. 
The K,CO3 remaining is 
0.58 + 1.53 = 2.11 mole 
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3.83 


Reaction (6) 
K,CO; dissociates at high temperatures: 


K,CO; = K,0O + CO, 
2.11 2.11 2.11 
From the above reactions, 7.27 mol No, 2.11 
mol K,0, 1.53 mol K.S, 0.77 mol S» and other 
gases CO, and CO are produced. In this case 
there may also be the water gas reaction, which 
produces H) gas. Therefore, in consideration of 
this reaction, we have: 


1 kg of mixture= xCO.+yH.+zCO+ 
u H,O + 3.64 No + 

0.77 S. + 2.11 K,0 + 

1.53 K,S + 0) 


From the number of atoms 


C=x+z=10.59 
O=2x+z+u+2.11 =23.31 
H = 2y + 2u = 3.00 


From these relations 


y=x-9.11 
z=10.59-x 
u=10.61—x 


(e) Calculation of Characteristic Values 


The reaction heat is obtained by the relation of 
the heat of formations: 


OQ, = Xco, + 29co + UF u,o) + 2.1 1dx,0) + 


1.5344.) —q+Aq 


From a rough calculation of explosion tem- 
perature, 7, = 2000 K was obtained. Therefore 
from Table 9 the equilibrium constant K(T) = 
4.46 was obtained. Therefore, 


z-u _ (10.59—x)(10.61-x) 


BES £(2 9.11) 


= 4.46 


From this equation we have 
x= 9.16. 

Further we obtained 
y= 9.16-9.11 =0.05 
z= 10.59 — 9.16 = 1.43 
u= 10.61 — 9.16 = 1.45 


We also have: 


Ix,0] = 86.8 kcal 
Ix,s] = 87.3 keal 
g = 7.27 dxnxo,) + 1-004 ¢1,0) + 9-1096.4,,05) 
= 7.27x119.54+1.00x 68.31+ 
0.10 x 226.8 
= 959.75 kcal 
Aq =0.57x(x+y+z+ut+3.64+ 0.77) 
9.16+0.05+1.43 + 
1.45 +3.64+0.77 
=9.41 kcal 


= 0st 


Therefore the value of Q; is obtained: 


Q; =9.16 x 97.80 + 1.43 x 29.70 + 
1.45 x 68.31 + 2.11 x 86.8 + 
1.53 x 87.3 — 959.75 + 9.41 

= 403.75 


From the internal energy of each product: 
Q, = xXEco, (1) + yEy, (1) + ZE eo (T) + 
UE, o(T) + 3.642, (1) + 0.77 Ey (1) + 
211E, o(T)+ 153.5) (T) 


[K,0] 


From of values of Q> are calculated from Ta- 
bles 6 and 8: [For £,, , the latent heat of evapo- 


ration (10.54 kcal/mol at 20 °C) is added. ] 


2.11E 


(T) 


[K20] 


A graph of T vs Q> is then drawn and the tem- 
perature at QO; = 404 kcal is read. This is the 
temperature of the explosion (7,). 


T, = 2050 K. 


The specific volume is: 
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Vo =(xt+ y+z+ut3.64+0.77)x 22.41 
oes 


1.45+ 3.64+0.77 
= 370 liter 


} 22.41 


The force of explosives f, is 


f, =0.3782V,T, = 0.3782 x 370 x 2050 
=0.287x10° dm 


4.3.3 Fuse, Shells and Pasting 


For the main fuse, industrial fuse of the 
second class was used. The shell was made of 
newspaper, which is the most rational choice 
because the breaking strength of newspaper is 
almost negligible. The internal diameter was 
105 mm, with a thickness of 2.5 mm. This is 
called four sun. Another shell called five sun 
had an internal diameter of 132 mm with a 
thickness of 3.0 mm. 


For the types of paper used in pasting the 
shells, the following symbols are used: 


Japanese paper Kraft paper 
Wie Ms 
Ws M, 


where W means Japanese paper and M means 
Kraft paper. The subscripts show the number 
of pasting layers of paper per a diameter length 
of sun (1.2 inches). For example, with the 
symbol W), for 5 inch shell (4 sun shell), the 
number of pasting layers should be 16 x 5 = 
80. Figure 14 (A) shows the layers of Ws and 
(B) shows the layers of Wi. for 5 sun (6 inch) 
shells. 


Figure 14. Sections of pasted shell. 


In this case sample shells were pasted with 
two-folded Japanese paper and unfolded Kraft 
paper. The total sum of the actual paper sheet 
area divided by the theoretical value S, was 
1.08 for Japanese paper and 1.17 for Kraft pa- 
per. 


The wheat paste used for preparing the two- 
folded Japanese paper was wheat starch and 
water in the ration 1:10. The paste used for 
pasting paper onto the shell was wheat starch 
and water in the ratio 1:3. 


The strength of the paper was measured ex- 
perimentally. The test pieces were made as for 
tests of metal tensile strength. The strip was 20 
mm wide. The sample consisted of no fold, 2 
folds, 5 folds and 10 folds of paper and the 
tensile strengths were measured. The humidity 
was 77% and the air temperature was 16.5 °C. 
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Table 10. Breaking Strength of Paper. 


Number of 
Paper Type Folds 


Japanese paper 
gundo 
M.C. = 8.8% 


Kraft paper 
(Cement bag) 
M.C. = 10.6% 


Newspaper 
M.C. = 10.9% 


= 
OON A AIN OO = 


M.C. = moisture content 


In Table 10 the symbol (—) indicates that the 
fibers are arranged at a right angle to the pull- 
ing direction, ( | ) indicates that the fibers are 
parallel to the pulling direction, and (+) indi- 
cates that the fibers are folded alternately (par- 
allel and at right angles). The tensile strength 
of (—) and (| ) are quite different from each 
other. The ratio of (—)/( | ) = 0.48 to 0.55. The 
mixed strength (+) almost equals the sum of 
(—) and (| ). The tensile strength of Japanese 
paper is almost the same as that of Kraft. They 
show 270-290 kg/cm’ in the direction of the 
fiber length and 140-150 kg/cm’ at right angle 
to the fiber. This is very small in comparison of 
the tensile strength of steel (3500-5000 
kg/cm’). The tensile strength of newspaper is 
25-35% of Kraft paper and very weak. The 
Japanese and Kraft papers have almost the 
same tensile strength. The thickness of the 
Japanese paper is 0.08 mm and the Kraft is 
0.14 mm. Therefore, two sheets of the Japa- 
nese paper correspond to one sheet of the 
Kraft. 


Fiber 
Direction 


Breaking | Tensile 
Thickness Strength 


(mm) 


4.3.4 Construction of the Sample Shell 


As sample shells, we used ring shells to 
make the analysis easy after the experiment. 
The construction of the ring shell is shown in 
Figure 15. The ring shell projects stars in the 
form of a ring. Cottonseeds were used to fill 
the space that would normally be occupied by 
other stars. 


A “Line of Stars” was drawn on the shell 
case to show the position of the ring. Three 
lugs were attached to the shell case to allow 
the shell to be suspended in a fixed orientation 
with respect to the plane of the star ring, as 
shown in Figure 20. 


The cottonseeds were put into a container, 
water was poured in to a known volume, and 
the volume of the cottonseeds was obtained. 
The true volume of 1 kg cottonseeds is 1.211 
liter. Therefore: 


oO 
A= 
V —ny, -1.211(a' +0") 


Page 36 Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


Photographing Direction 
Cotton Seed 


Paper Bag 
Shell ; 
Prime 


Star 


Burst 
Charge 
Black 


Tube Powder 


(A) Sectional View 
(horizontal) 


Figure 15. The construction of the sample ring shell. 


where 


the quantity of the burst charge 
without the binding material (kg) 


v, | the volume of a star (I) 


the weight of the cottonseeds for 
the burst charge (kg) 


, | the weight of the cottonseeds for 
packing material (kg) 


A | the loading density (kg/l) 


The calculated data are shown in Table 11. 


4.3.5 Table of Sample Shells 


The sample stars manufactured by the above 
method are arranged in Table 11. The symbols 
used in the table are as follows. For example, 
D,BsHWs means 


Line of 
Stars 


(B) Sectional View 
(vertical) 


4 sun shell (5 inch shell) 
8 mm blue star 
Burst charge: haisan (potassium 


chlorate burst charge) 
Japanese paper sunto 8x4 folds 
pasting. 


Shells 1 to 48, which all contained blue 
stars, were used to determine the effect of the 
movement on stars of: the size of the shell, the 
size of the stars, the type of burst charge, the 
kind of paper pasted, and the number of pasted 
layers of paper. 


Shells 49 to 56 were used to determine the 
effect of velocity on the star. Specifically, 
shells 49 to 53 contained three types of stars in 
the same shell, and were used to compare the 
velocity of stars without the effect of manufac- 
turing deviations. Shells 57 to 60 were so called 
manboshi (ordinary chrysanthemum), in which 
the shells were filled with stars. This was for 
comparison with the ring shells. Shells 61 and 
62 were commercial shells that were compared 
with the test sample shells. 
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Table 11. Sample Shells. 


Burst Charge 


Shell No. Quant. | Net Wt. | Cottonseeds 
Dia Color | (mm) | Stars | Form | Type (kg) (kg) B. Chg. (kg) 


i00 | G0 $00 }0@0 00 |00 00 }~@ ia 


P 
H 
Ss 
P 


VVVVUA WH H|= 
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Table 11. Sample Shells. (Continued) 


Cottonseed Loading Density Paper Pasting 


Filling of Burst Charge 
(kg) (kg/l) 


| DaBr2PWie | 
| D.Bi2PM, | 
Ds B 12 P M 8 
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Table 11. Sample Shells. (Continued) 


Burst Charge 


Shell Dia. No. Quant. | Net Wt. | Cottonseeds 
Dia. | Color | (mm) | Stars | Form = (kg) (kg) B. Chg. (kg) 


wf DsBr2SWe tL | blue | 12 | 34 | Ring | 

a DsBi2SWie Ld | blue | 12031] Ring 

[DsBieSM4 Lf blue | 120) 34 | Ring | 
D;Bi25M, 

wf DsBr2PWe tL | blue | 120] 34 | Ring | 

5. DsBi2PWie Ld blue | 12031] Ring | 

[DsBi2PMg Lf blue | 120) 34 | Ring | 

D 5 B 12 P M 8 

wf. DsCi2G12Bi2SWe fd | 12 | 30. | Ring | 

) pDsG12SWe 12 30 | Ring | 

Ose Wg dl blond nt | Ring | 

D;C,.SW, 

jDsCi2Gi2Bi2PWe | tt | 12 | 30 | Ring | 

fj DsGi2PWe tt 1230 Ring | 

DeGeP Wee a 8 | Ring 

D 5C 12 PW; 

of DSBS Wg cl ed Ef 860, | Full | 

3 DsG12SWe 12 318 |Full 

[DsBi2SWe | Asun | 120 | 190 | Full | 

D.Bi25W, 

D,B,.SW, 

D,B..SW, 


Note: Shells 57-60 were chrysanthemum shells, and Shells 61 and 62 were commercial shells. 


ViV|DD:'D:'D/ VV VV: WW! VV VV WWW 


— means that the color of the stars were not clear. 
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Table 11. Sample Shells. (Continued) 


Cottonseed | Loading Density of Paper Pasting 


Filling Burst Charge Type No. of | da-di Total 
(kg) (kg/I) of Paper | folds | (mm) | Wt. (kg) 


wf DsBr2S Wed. 

, DsBi2SWie dn. 

[OSB SMa da 
D;B,.SM, 


[DsCi2Gi2B12PWe | 
D;G 12 PW; 


Note: color of stars: 


Shell 49 Amber to Green to Blue CGB 


Shell 62 Double petalled xX 


Shells 61 and 62 were commercial shells. 
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Figure 16. Revolving shutter. 


4.4 Preparation for the Test 


4.4.1 Main Instruments 


(1) Camera and Accessories: Camera (focus: 


180 + 1 mm, f= 4.5) with cabinet-type 
dry plate. 

(2) Revolving Shutter with a Synchronizing 
Motor: 25 revolutions per second, one 


wing with shutter angle of 60° (Fig- 
ure 16). 


(3) 16-mm Movie Camera: 128 frames per 
second, | inch focus, f= 1.8, 16-mm Ko- 
dak film (100 feet) (Figure 17). The re- 
cord from this camera was used only for 
reference. 


(4) Iluminators 

(5) Surveying Instruments 

(6) A Recorder 

(7) Distance Poles with Lights 


(8) Electric Wire for Revolving Shutter 
(9) Side Camera (not special) 


140 
130 
120 
110 
100 


Frames per Second 


K<e) 
oO 


Time from Start (s) 


Figure 17. Number of frames per second as 
the time passes. 
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4.4.2 The Arrangement of the Instruments 
for Burst Charge Experiments 


The installation of the test shell is shown in 
Figures 19 and 20. The plane of the star ring, 
which was marked on the outside of the test 
shell (Figure 15, (B)) was placed vertically so 


> e A 


a 16°6' 30" 
Position of 
the Test Shell 


(H=0) 


Side Camera 


No. 1 to 7: Light Pole 


that it coincided with the vertical plane includ- 
ing the distance poles with lights (Figure 18). 
The test shell was held firmly with three pieces 
of wire 1 meter above the ground (Figures 19 
and 20). 


Revolving 
Shutter 
+ Lens 
= 
oO 
8 
Oo 
© 76°15' 30" Main Camera with « Shell” 
y Revolving Shutter P eee of 
SV Auxillary explosion) 
< Explosion 


Mark (H=0) 


11°46'30" 


/ 


35 mm Movie Camera (H=19.75 m) 


Figure 18. Arrangement of the instruments for the exploding shells. 


Figure 19. Installation of the test shell (1 of 2). 
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ae ws ‘S Si. eas: 
Figure 20 Installation of the test shell (2 of 2). 


The distance poles with light were set at in- 
tervals of 20 meters from the bursting point 
toward the west (No. 1-6). Another distance 
pole (No. 7) was set at 60 m from the explod- 
ing point toward the east. Two other distance 
poles were set on either side of the explosion 
site to the north and south to mark the horizon- 
tal line for the side camera. All of the lights on 
the poles were installed so that they were level 
with the test shell, which was hung | meter 
above the ground. The distances between the 
poles (small electric bulbs) were measured 
with a tape measure that was verified with a 2 
m steel tape measure. 


The height of the lens of the revolving shut- 
ter camera was a little higher (0.64 m) than that 
of bursting shell. However it was regarded as if 
they were on a same level. A transit was placed 
at the No. 2 pole and the axis of the lens of the 
revolving shutter camera was adjusted so that 
the axis fell at a right angle to the pole line. 
Thus these installations could catch the image 
of about 1/4 of the chrysanthemum. 


Small electric bulbs were used for lights on 
the distance poles for the preliminary test, but 
from then on candle lights were used to save 


the electric cells. Only the light on pole No. 2 
used both an electric bulb and a candle. How- 
ever, based on the results, only the candle light 
was useful. 


4.5 Progress of the Experiment 


The first day (6 November 1954) 


Preliminary tests were conducted. The test 
shells were burst on the ground, the bursts 
were photographed, and the film was devel- 
oped to confirm the results. Time: 9:30 to 
10:00 PM, 0 °C, 100% humidity. 


The second day (7 November 1954) 


Regular tests were conducted. Ignitions with 
electric devices did not work well and failures 
occurred with abnormal delay. Then we used 
black match with good results. Time: 7:15 to 
10:15 PM, 80% humidity. The voltage was too 
low to power the sound recorder. 


The third day (8th December 1954) 


Time: 6:30 to 11:40 PM, Windy 3—5 m/s, 
light rain. 


Operations 


The shells bursts were photographed by the 
camera with the revolving shutter and by the 
16-mm movie camera. For confirmation, the 
side camera was used to show whether or not 
the stars filled a vertical plane. 


The tests were conducted as follows: 
e The second day: Shells 1-28 


e The third day: Shells 29-60, also Shells 2 
and 8, which were left over from the sec- 
ond day. 


Photographs for the following tests could 
not be obtained: Shells 1, 12, 46, and 52. 


The photographic results from the 16 mm 
movie camera were not easy to use because of 
the deviation of delay times of the ignition of 
the test shells. However, the following shells 
were caught almost perfectly by the camera: 
22, 28, 33, 34, 36, 37, 40, 42, 43, 44, 45, 47, 
49, 52, and 53. These photographs do not have 
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So 0 


Figure 21 A. An image taken through the 
revolving shutter camera from the front. 


enough accuracy to determine the velocity of 
the stars, but they are useful to explain the way 
the stars move near the explosion and the way 
the burst charge burned. 


The cycles per second of the electric power 
that was used for the synchronized motor of 
the camera, from the report of Itsukaichi Power 
Plant, were as follows: 


| December_|8PM__| 11PM | 15 PM_ 


The following photos show examples of the 
results obtained from the above experiments. 
Figure 21 A shows an image taken by the cam- 
era with the revolving synchronic shutter. The 
locus of the stars appears as a series of dotted 
lines. Figure 21 B shows an image from the 
side camera, Figure 21 C shows images from 
the movie camera. All of the photos were taken 
through an open iris. 


Figure 21 C. The images taken by the 16-mm 
movie camera. 


Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. Page 45 


Chapter 5 — Supplementary Experiments: Measurement of 
the Burn Velocity of Burst Charges and Stars 
in the Normal Atmosphere 


To obtain these data, the 16-mm movie camera was used. A Black Powder pasted paper strip was 
arranged on a piece of metal mesh. Grains of burst charge, or stars, were arranged on the strip at inter- 
vals of about 5 cm, and the paper strip was ignited at one end. The grains or stars ignited one by one. 
They were photographed from a distance of about 1.5 meters. 


Table 12. Burn Velocity of Stars (Measured by 16-mm Movie Camera). 


a aa Number of Frames 
Colored sper 


Note: The figures in () show the average burn time when the speed of the camera was 120 frames per second. 
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Table 13. Burn Velocity of Burst Charges. 


H 


P 


Ss 


No. Frames 


No. 


No. Frames 


No. 


No. Frames 


Note. The figures in () show the average burn time when the film speed was 120 frames per second. The photo- 
graph did not show clear images of the starting points and ending points. The accuracy of the data is therefore 


not very good. 
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Chapter 6 — Formulae for Calculating the 
Velocity of Stars 


First, the author studied the law of motion 
of stars in the air. The velocity of stars at dis- 
tant points from the shell burst was then re- 
lated to the initial velocity of each star. The 
conditions that produced these initial veloci- 
ties were investigated, and from this informa- 
tion, equations useful in the designing of 
shells were developed. 


6.1. Arrangement of Experimental 
Data 


6.1.1. Photography with a Revolving 
Shutter Camera 


The scale of the image taken on the dry 
plate of acabinet was 1/757.6. This is based on 
the distance between the calibration poles. The 
images of the star flight lines were numbered 
from the lowest trajectory. For example, Shell 


: 
Lal 
ae oo te Sess Sn 


cooreee 


No. 231 represents the lowest trajectory of the 
star from the Shell No. 23. 


6.1.2. Photography with a 16-mm Movie 
Camera 


The calibration poles were not clearly ob- 
served with the 16-mm movie camera, how- 
ever, shells Nos. 22, 37, 49 and 51 barely in- 
dicated the position of the light of the calibra- 
tion poles. The scale of the image was 1/142.9. 
Figure 22 shows the image from the 16-mm 
movie camera. 


6.2. Trajectory of the Star in the Air 


From the trajectories in the photographs, it 
became clear that the velocity of a star is 10— 
100 m/s, which is a fairly low velocity com- 
pared to the speed of sound. The simple ap- 
proach to calculating air resistance outlined in 


Figure 22. The trajectories photographed by the 16-mm movie camera. 
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Chapter 3 should therefore be adequate. In 
other words, it should be possible to use equa- 
tion 38 for moving stars. To solve this equa- 
tion for moving stars, it is sufficient to con- 
sider only the horizontal trajectories for vari- 
ous examples. 


The above equation was modified to a 
more general formula: 

2 

xk y 0 
dt pt) 

From the results of the experiments, the 
values of n and k are obtained as follows. The 
sample stars used for these calculations should 
be selected so that the star’s flight is as hori- 
zontal as possible and has a minimum devia- 
tion from the line of calibration poles. Flying 
stars are regarded as being free from the ex- 
plosion gas. 


6.2.1. The Reaction of Air upon the Star 


The following calculation is for Star No. 
481. 


As described earlier, the trajectory of a star 
appeared on the photograph as a dotted line. 
Table 13 shows the order of calculation along 
the dotted line. The number on the left side 
shows each dot from the start on a trajectory 
that had the lowest horizontal flight. Number 
1 means the dot nearest the burst point, 2, 3, 4, 
etc. are the order of the star’s position as it 
traveled out from the center. The value of x is 
a coordinate on the x-axis associated with each 
dot numbered from the burst point. The value 
of Ax is the distance between | and 3, 2 and 4, 
3 and 5, etc. (See Figure 23.) 


Table 13. An Example of Calculations for a Star That Flew Horizontally. 
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These were recorded at the heads of each in- 


Dots a. Se termediate dot (Figure 23). The time interval 

| <= aK | between each pair of dots is 2/25 seconds. 

| < Ax | Therefore the value of Ax divided by 2/25 s de- 

| < AX > | notes the negative acceleration. The symbols v 

Figure 23. Diagram of variables in Table 13. and dv/dt are used in place of 7 and aw in 

For example from Table 13: the following paragraphs to avoid any confu- 

21.84 — 15.96 =5.88 sion. Table 14 shows the values of such data 
14.32 = 18.92 _ 5 40 obtained from the test shells. 


26.76 — 21.84 = 4.92 


Table 14a. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14b. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14c. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14d. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14d. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (Continued) 


Table 14e. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14e. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 


8 


©) 10) : 0) }O) $00 00:0 
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Table 14f. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14f. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 


Table 14g. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14g. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 
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Table 14h. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14i. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14j. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 


58.8 — : 54.6 : 65.0 
53.9 119 : 51.6 ; 60.1 
49.3 106 : 49.0 : 56.0 
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Table 14j. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 


Table 14k. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 


—_ 


© }00 
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Table 14k. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 


Table 141. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14]. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 
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Table 14m. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14m. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 
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Table 14n. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 14n. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 
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Table 140. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. 
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Table 140. Horizontal Velocities and Horizontal Negative Accelerations of Stars that Flew 
Almost Horizontally. (continued) 


6.2.2 The Burn Time of Stars 


The burn time of a star, not including the 
core, and initial ignition material is determined 
by counting the dots along the trajectory in the 
photograph. The number of dots is shown in 
Table 15. 
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I. The Stars of Bg. ** n denotes the number of dots counted along 
Table 15a. Number of Dots in the Trajectory the trajectory in the photograph. 
of the Stars. ea The first 1 or 2 digits (e.g., 28) denotes the 
[No* | n* | shell number; the second | or 2 digits 
Ne pe fo Nef (e.g., 10) denotes which star this is from 
the horizontal position. 


II. The Stars of By. 


Table 15b. Number of Dots in the Trajec- 
tory of the Stars. 


No. denotes the number of the star. 


** n denotes the number of dots counted along 
the trajectory in the photograph. 


= No. denotes the number of the star. 
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Table 16. The Burn Time, Burn Velocity and the Probability Deviation of the Burn Time (from 


the Image of the Trajectory of the Star). 


ane te A | le aL eee 


Ill. The Stars of G12 


Table 15c. Number of Dots in the Trajectory 


of the Stars. 


* No. denotes the number of the star. 


we denotes the number of dots 
counted along the trajectory in the photograph. 


The image of the trajectory of the stars in 
the photograph shows that at the start, the light 
output from the burning of the Black Powder 
type ignition composition was weak (Photo 8A). 
The revolving shutter did not capture the image 
of the trajectory at the start of the ignition. 
However, at the start of the color composition, 
the image was captured sharply. At the end, a 
core flash appeared. In this case, the light is too 
bright and could not be captured clearly. Con- 
sequently, the burn rate of the star can be calcu- 
lated only during the color part of the composi- 
tion. 


Therefore, for the calculation of the burn 
time of the color part of the composition of a 
star, we can use the number of dots in Tables 15 


a, b, and c. The burn time of the color part of 
the composition is therefore 

xn 1 _ 

T, =—x— =0.04n 

N 25 
where 7, is the average burn time for N stars 
with the color composition, 1 is the number of 
dots, and 7 is the average number of dots per 
one star. The probability deviation of the burn 
time is 


Snap 
7 = 0.6745 x 0.04, |; —=—____ 
N-1 


The burn velocity w; is as follows, when the 
burning proceeds from the surface towards the 
core 


where R; is the radius after the ignition layer 
has burned off, and FR, is the radius of the core. 
Table 16 shows the results of the above calcula- 
tions. 


From the above descriptions, the burn data 
of stars were obtained by means of three meth- 
ods: stopwatch timing of stationary stars (4.3.1), 
from the 16-mm movie camera with stationary 
stars (4.3.1) and from the revolving shutter 
camera for moving stars. They are summarized 
in Table 17. 


Table 17. Comparison of the Data for Burning Stars. 


Stationary 


Stopwatch (s) 


16-mm Movie Camera (s) 


Revolving Shutter Camera (s) 
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The data from the stopwatch shows the total 
burn time minus the burn time of the core. The 
data for the burn time of stars are not very dif- 
ferent between stationary and moving stars. 


6.2.3 The Sectional Density of Stars 


Since stars are burning as they fly through 
the air, the sectional density changes with time. 
The sectional density of the star at the nth dot is 
calculated as: 

2 = 
where m, is the mass of the star at the nth dot 
and o, 1s the sectional area at the dot. For con- 
venience, the number of the dot is applied as v 
= 0, 1, 2, 3, ..., 18, 19, 20, etc. in the opposite 
order as the dot number 7 = 20, 19, 18, 17, ..., 
3, 2, 1,0, 


etc. When the burn rate is kept constant, the 
sectional density p/ is denoted as 


R,+vw ; 3 4 
+4( 4 45) Ri} tap 
0, uaa . 
25)" 


where m1» is the mass of the core. 


For blue stars (Bg and B;2) and green stars 
(Gi), the values of p’ were calculated as they 


are in Tables 18 and 19. 
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Table 18. Sectional Densities of Stars, g/cm’. 


I. Blue Stars (Bg and B;2); w,/25 = 0.078 cm/s, 
5 = 1.62, Ra= 0.25 em 


II. Green Stars (G2), w,/25 = 0.047 cm/s, 
6 = 1.64, Ro = 0.25 cm 
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6.2.4 Calculation of a Function, F(v) = kv". The value of —(dv/dt), the value —dv/dt in 
Table 14 are usable. For the values, the number 
v should be used. For example, Star No. 481 is 
used as the following example. 


In the equation of motion 40, the air resis- 
tance of kv" is calculated by putting (dv/df), to 


dx and p’ to ‘(t) 
dt? aa 
dv 
—| — |p =k" =F 4] 
Ga V(Y) (41) 


— 
oe 
ee 
ee 

e 

e 

e 
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The notation 1 indicates that the characteristic of the logarithmic value is —1. 


For example, if the number Z = X x 10”, where X is > 1 and <10, then the log (Z) is written as “Y.(logX)” 
where the dot is a decimal point. The characteristic (the number to the left of the decimal point) shows the power 
of 10 (i.e., Y), while the mantissa (the number to the right of the decimal point) shows the power to which ten 
must be raised to give X. If the number Z is less than one, the characteristic will be negative and is written with 
a bar over it. The mantissa is always positive. 


The values v, were taken from Table 14. The 
relation of the stars concerning log v, and log 
F(v) are shown in Table 19. 
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Table 19a. The values of log v, and log F,(v). log FQ) log F,(v) 
v_ | logv, 100 logv,, 100 


| | Shell No. 51 Shell No. 62 


BPP 


Shell No31 | No. 31 Shell No.41 No. | Shell No.41 | 


Table 19c. The values of log v, and log F\(v). 


eas, i 7 on : 7 
log log 
v_ | logv, 100 logv,, 100 


| | Shell No. 72 Shell No. 82 


pala 
1 
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Table 19d. The values of log v, and log F,(v). Table 19e. The values of log v, and log F\(v). 


| i, (v) 
oO ae, ES a 
ea tt) Sa at) v_ | logv, 2 100 logv,, 100 


| ShellNo.91 | ~~ ShellNo. 101 | | | Shell No 112 Shell No. 132 


2 
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Table 19f. The values of log v, and log F,(v). Table 19g. The Values of log v, and log F,(v). 


| F(v) 
og —-—— 
See Lt) SSE a) v_ | logv, : 100 logv,, 100 


| Shell No. 142. | ShellNo. 153 | | | Shell No. 162 Shell No. 172 


(j=) 


10 [NIi@ i iO lo iN! 


jo) 


oe. 
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: ar 
eel 
ef 
eel 
“eae 


© 100 100 100 100 10 IN SENE: 
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Table 19h. The values of log v, and log F,(v). Table 19i. The values of log v, and log F\(v). 


FQ) FQ) 
log —— log 
See tt) Sa at) v_ | logv, 100 logv,, 100 


| Shell No. 182. | — ShellNo.191 | | | Shell No. 202 Shell No. 212 


POT 7O1 iW iW: 


Ceacoe | 
ING: 


ae 
2 
Pe 
He 
et 
et 
et 


© & 1 1 > ON NIN: 
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Table 19j. The Values of log v, and log F,(»). Table 19k. The values of log v, and log F,(v). 


FQ) FQ) 
log ——— log 
Sei tt) See a) v_ | logv, 100 logv,, 100 


| Shell No. 222, | Shell No. 232 | | | Shell No. 241 Shell No. 252 


(j=) 


alo! wis 


oa 


o 


Heo ior) 


(=) 
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Table 191. The values of log v, and log F,(). tog (”) gow 
v_| logv, 100 logv, 100 
F F 
log— () log — (¥) 
v_ | logv, 100 logv,, 100 


| | Shell No. 261 Shell No. 272 


| | Shell No. 281 Shell No. 291 


F.(v) 
log 
v_ | logv, 100 logv,, 100 


| | Shell No. 302 Shell No. 312 


Table 19m. The values of log v, and log F,(). 


el =e 
log log 
v_ | logv, 100 logv,, 100 


| | Shell No. 281 Shell No. 291 
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Table 190. The values of log v, and log F,(v). Table 19p. The values of log v, and log F,(v). 


FQ) FQ) 
log ——— log 
Sone itt) See a) v_ | logv, 100 logv,, 100 


| Shell No. 343 | Shell No. 352 | | | Shell No. 362 Shell No. 371 
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Table 19q. The values of log v, and log F,(v). Table 19r. The values of log v, and log F,\(v). 


| i, (v) 
Oo os 
Snes Lo) See at) v_ | logv, 100 logv,, 100 


| Shell No. 382. | — Shell No. 392 | | | Shell No. 403 Shell No. 411 


oa 


oO 
POPNE OO: AR WENI>I0 


7 
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Table 19s. The values of log v, and log F,(v). Table 19t. The values of log v, and log F,(v). 


FQ) FQ) 
log ——— log 
= tt) logv,, a) v_ | logv, 100 logv,, 100 


| ShellNo.422 [~~ ShellNo.442 | | | Shell No. 452 Shell No. 471 
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Table 19u. The values of log v,, and log F\(v). log FQ) log F(v) 
v_ | logv, 100 logv,, 100 
fi i | 


Shell No 487] No. 481 Shell No. 502] No. 502 


| | Shell No. 481 Shell No. 502 


Bt v_| logv. 100 logvs. 100 


ae 
181 | |__ Shell No. 503 Shell No. 505 
8 Na ee aS ee SE 


i) 
K 
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FQ) 
log log 
v_ | logv, 100 logv,, 100 


| | Shell No. 503 Shell No. 505 
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6.2.5 Calculation of the Values nm and k for 
the Motion Equation 40 
SNe Sa I ae Sco on 4 


| ShellNo.542. | ShellNo.543_—i| A graph of the relationship of log v and log 
F‘(v) may be obtained as in Figure 24. 


'<— log F(v)-2 


log F(v) 


a 
log v-logv, 


log v 


Figure 24. A graph for the relationship 
between v and F(v). 


As a trial, a line AB is set as: 
log F(v)-2 
logy —log v, 


where n' is the slope of AB and v, is the value 
at log F(v). From the above equation we have: 


log F(v) =n'logv+(2—n'logy,) (42) 
On the other hand from equation 41 
log F,,(v) =nlogy,, + logk (43) 
Comparing equation 42 with 43, the next re- 
lations are obtained, 
ni=n 
(44) 
2—-n'logy, =logk 


The values n' and log v, can be obtained by the 
method shown in Figure 25, and the values of 
and & can be calculated. Figure 25 shows the 
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relations of log v and log F(v) = log(F,(v)/100) 
+ 2 (Table 19). 


3.2 T T T T T 
r © 051 i 
307 @ 222 a 
r p:] 
96 © 362 a J 
| A 403 age 
26 [ & 542 og | 
_ 244 og 7 
= PP iga 
= 2.2 i Sie 4 
ae) A 
2.0 b oes J 
rs 
1.8 + S at 4 
r é 
1.6 + ul? 4 
L ae 
144 # J 
1.2 I a, | l L i l | 1 
08 10 12 14 #16 £18 
log v 


Figure 25. The relationship of log v and 
log F,,(v) (See Table 19). 


The values n’ and log v, and their average 
values and the probability deviations are shown 
in Table 20. 


Table 20. The values 7’ and log v,. 
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In the Figure 26 the values of n for Bg are 
biased toward lower values, and B;2 is biased 
toward higher values. However, it could be as- 
sumed that they fall in line with each other. 


1.4 T T T T T T 
1.34 ; ae | 
° fo) 
o 12> Be = . Bey e | 
> o *° $ e 
D 1.1 4 
2 @ B8D-4 
1.04 O B-12 D-4 
@ B-8D-5 
0.9 + © B-12D-5 
@® G-12D-5 
0.8 


a 
1.6 18 20 22 2.4 
n' 


Figure 26. Ballistic constants log v, and n'. 


6.2.6 The Influence of Weather Conditions 
on the Ballistic Constants 


On the first day of the test, it was calm with 
no wind. On the second day, it was windy (4—5 
m/s) and a little rainy. If there were large influ- 
ences on the values of n and log v,, they would 
be divided into two groups. Figure 27 shows 
the values for n and log v, obtained on the first 
and second days using different symbols. The 
values of n and log v, of the first day appeared 


1.40 
1.30 

> 120 Dube 

D- . 

2 o's . 
1.10 be 


to be lower than those of the second day. How- 
ever, for simplicity, these differences were ig- 
nored because such influences would probably 
not occur in practice. Therefore the average 
values of n and log v, were used for the follow- 
ing calculations. 


6.2.7. The Differential Equation Describing 
the Movement of Stars 


From equation 44, the values of ” and k are: 


n=n'=1.93 
log k=2—n' log v, = 2 — 1.93x1.20 = 1 .684 


Publisher’s note — there was no equation (45) 
From above 
k = 0.483 


The differential equation for the horizontal mo- 
tion of a star in the air should therefore be written 
as: 


d’x 
dt? 


p(t) = = 0.48308 (46) 


6.2.8 Formula for a Star Moving in a 
Horizontal Direction 


From equation 46: 


ee (47) 
dt pt) 
3 Pol oo 


e : First Day (1954, 12, 7) 
© : Second Day (1954, 12, 8) 


AS) SO. oe 8 


gs Reese 0 | 


2.2 2.0 24 2.5 


n 


Figure 27. Influence of the weather on the ballistic constants. 
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where k = 0.483 and n = 1.93 


Integrating equation 47 we have: 


Se = (48) 


pS t at 
nv Sop) 


Here, v, is the velocity at some start time 
(t= 0), and we can select any convenient value. 
Therefore, v, is taken either as the initial veloc- 
ity of the star or the velocity at the color change 
of the star. Solving equation 48 we have: 


raft n)(-k) ['—@ c (49) 


Introducing 


_ ax 
dt 


we have: 


y 


= ffm" +(1—n)(-k) a “dt (50) 
0 oD 


Equations 49 and 50 are the fundamental equa- 
tions to express the relations tv and +x for 
horizontal movement. 


(1) For Stars That Consist of a Single 
Composition 
Using the symbols: 


t time 

m the mass at time ¢ 

r the radius at time ¢ 

R the initial radius at time t = 0 
p(t) sectional density at time f 


we have: 
m 
p(h=—; (51) 
amr 
On the other hand, 
m= eee 6 
3 
where dis the density of the star. 


Therefore, 


r') =515 (52) 
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Introducing the following formulas: 


r=R-oat 

R (53) 
o=— 

T 


where T is the total burning time of the star, we 
have: 


y=a(I-2}=R(1-2) (54) 
where 

t ' 

Cages (54') 


The symbol 7 denotes the ratio of the burn time 
t to the total burn time T. 


Next, the following calculation is carried out: 


eat dt 23 Doe dt (55) 
opt) 4Rd-° it 
i 
Introducing equation 54' and 
1 
dt=—dt 55' 
7 (55') 


into equation 55, we have: 
he 3 T ¢*_dt 
4 Rd*° (1-7) 


3 T 
=——-2.303—log(l- 56 
i Plas. (56) 


Introducing equation 56 into equation 49 we 
have: 


a. 
l-n 


vy +(1 =n)(-#)x{ -332.308 
v= 
T 
= log! —7) 
RO a 
ry (1=n)(k)x(-8%2.308) ra 
4y," 


=V 
Ss log(1-7) 


Introducing a new symbol into equation 57 
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(1—n)(—k)x(-3x2.303) 7 


A,= = 
4y,” Ro 


AY 


(58) 


and calculating it with n = 1.93 and k = 0.483, 
we obtain: 


0.93 
A =0.774% + 
| Re (58') 
= 0.7740 1. 4 
RIT 6 


The value of A, is almost proportional to the 
value of v, and inversely proportional to the burn 
rate R/T and the density 6 of the star. 
From equation 57: 

1.075 


v=v[ (I-A, log(1-r)} See 


For comparison with the results of the follow- 
ing section, it is convenient to write this as: 


vw) 41-($ Jios(t-2)' dr (59) 


Considering equations 56 and 59 we have from 
equation 50: 


x= wr, 4! -(4 tos} dr (60) 


Equations 59 and 60 are practical equations to 
calculate the motion of a star that consists of a 
single composition. 


(2) For Stars That Consist of More Than One 
Type of Composition (Color Changing 
Stars) 

In this case the layer of each color has an in- 
dividual burn time (7;), outside radius (R;), and 

density (6;). (See Figure 28.) 
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Figure 28. A color changing star. 


First, the value of the sectional density is 
denoted as: 


pi)=—3 (61) 


The value of m is denoted as: 
4 (7° -R)d,+(R-R;) 6. 


2 
m= -1 (62) 
3 |+(R}-R;)6,+ ... 


—-R 
r=R -@t=R,- =f 
T, 
(63) 
afi 
Ls (R, /R,) 
Setting an equation as: 
T, 
T, =—___ (64) 
'  1-(R,/R,) 


Equation 63 becomes 


t \ 
=R|1-— |=R, (1-7 65 
( | (1-7) (65) 
here t, =¢/T,;. When R> = 0, equation 65 coin- 
cides with equation 54, and it becomes the case 
of the single colored star, where 7, =T,. In the 
same manner as shown in the previous section, 


the following integration was carried out. 


t tar’ 


Prowl 
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When the value m (see equation 62) is intro- 
duced into the above formula, we have 


r=f! mr°dt 
°4 [(° -R)6,+(R} -R3)6, + 
SHE. 
3 |(R} -R})6, +... 


When the next symbol c is introduced into 
the above equation: 


(6, -5,)R3 + (5, -5,) B+ a 
~|Op8 RG Fic 


we have 
2 
r=2f is dt 
4° 0rd, +¢ 


Equation 65 is introduced into the above equa- 
tion 


: si * TR? (1-2') dr 
aa" R(1-r) +e 


This integral is evaluated by making the substi- 
tution 


y=(l-ry 


dts -3(l —t)°dy 


which leads to 


j TR? (1-z') dt 
R(I-r) d+¢ 
“1 TRdy 


3R iyte 
oe 
3° "4 RSy+e 


"p2 
=-545 In(R°d,y+c)+C 


3 
ilies | 


= -> 52.303 og(( —r)+b)+C'! 


is 


where C and C’ are constants of integration, and 
b is given by 
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(67) 
The definite integral from r= 0 to r= T is 
then 
2.303 7, (l-ry +b 


T= 
4 Ro, 1+b (68) 


If the density is the same for all layers, 
6, =0, =0, =... 
then the value of b becomes 0 (b = 0). When the 
density of the star increases toward the center: 
0; 2 OO, Sas 
then we have b > 0. 


On the contrary, when the density decreases 
towards the center: 


0, > 0, >0,>... 
then we have b < 0. 


When the density does not change through 
all the layers, the value of b = 0, and we have 
the same equation as equation 65. When equa- 
tion 68 is introduced into equation 49, we have: 

1 


2) len 
x 


ve" 4(1=n)( x : 


, 3 
TF gg +b 
Ro, 1+b 


v= 


, (69) 
‘a (1=”)(-k)x(-2.303) | Te 
Avi," 


T, lo (i-z) +b 
Ro, 1+b 


=V 


In equation 58, a symbol A, was set as: 
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fe Meare), Fr (70) 
Ay,” Ro, 


Introducing n = 1.93 and k = 0.483: 
0.93 1 1 ' 
A, =0.774 vy? —_-— (70') 
R Eg 6, 


Then from equation 69: 


—1.075 


v=vy41 (4)oe 2) i (71) 


1+b 


When b = 0, this equation coincides with 
equation 59. And from equation 50, considering 
equations 68 and 70, we have: 


(1 i b —1.075 
_ et A, = te. + ; 
e=VE |. f c }toe ioe dt, 


(72) 


This equation perfectly coincides with equa- 
tion 60 when 5 = 0. 


Equations 71 and 72 are practical equations 
to obtain the relationship of v — ¢ and x — t when 
the first layer of the star is burning. To obtain 
the relationship of v — ¢ and x — ¢, when the sec- 
ond and subsequent layers are burning, first 
calculate the final velocity of the first layer. 
This velocity is then regarded as the initial ve- 
locity of the second layer v, and in place of 
R,, T,, and 6,, use R,, T,, 5, and equations 71 
and 72 to calculate the values of v and x. This 
repeats for the next layer, etc. 


6.2.9 The Curvature of the Trajectory of the 
Star 


The curvature of the trajectory of a star 
should be reasonable for our purpose. Therefore, 
when we design a chrysanthemum, we should 
obtain a proper value of the radius of curvature. 


Symbols: 
p radius of curvature 


mass of star 


m 
g acceleration due to gravity 
v velocity of star 

) 


inclination of trajectory 
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The centrifugal force must be equal to the 
centripetal force. Therefore, 


2 
mv 


=-—mg cos0 


From this relation: 


2 
pcos =—~— (72') 
& 


This equation gives the radius of curvature 
at any point of the trajectory. (See Figure 29.) 


Figure 29. The curvature of the trajectory of a 
Star. 
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6.3 Initial Velocity of Stars 


The cutting of the image of a flying star by 
the revolving shutter is not good at the start be- 
cause of interference from the sparks of the 
prime, which covered the surface of the star. 
Therefore, the initial velocity of stars should be 
calculated from where the image of the stars 
first becomes clear, near the bursting point. 
(See Figure 30.) 


6.3.1 The Formula of Initial Velocity 
Deduced from a Dot near the Bursting 
Point 


The equation of the motion of a star can be 
written as: 
d : kv 
ue gsin@ (73) 
dt p' 


dv . Rul : . 
where ei is the acceleration in the direction of 
t 


the movement, g is the acceleration due to grav- 
ity, A is a constant, and p' is the sectional den- 
sity of the star. To solve equation 73, the fol- 
lowing assumptions could be made because the 
initial calculation points are so near the bursting 
point: 


1) the acceleration due to gravity only affects 
the trajectory in the vertical direction, not in 
the horizontal direction. 


Bursting 
Point 


2) the constant k does not change during the 
motion. 


3) the sectional density does not change during 
the motion. 


With these assumptions, the equation of mo- 
tion can approximately be divided into two 
parts: 


the horizontal direction: 


Seen y (74-1) 


Tecpel yy (74-2) 


Author’s note: 


In this case the author could not find a good 
dynamic calculation for the motion of the stars 
because the trajectories of the stars near the 
bursting point are invisible, as they are in Fig- 
ure 30. Therefore, the following is based on an 
approximate calculation. 


From equation 74-1 we have: 


IE Invisible Region 


Figure 30. The image of stars near the blasting point on the photograph. 


Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


Page 93 


Eom aeaa Or: Table 21. Velocities of Stars near the Burst- 


dv a a a ing Point. 


Vv. P 


Integrating this equation: 


P (75) 


Integrating equation 75 we have: 


2 2 ry & ry 
Kat 2 vse 2 el (76) 


The initial velocity of the star can be written 
approximately as: 


Ve=Vi+Vy 


Therefore, the combination of equations 75 and 
76, gives us: 


oF y aky gp' oF 
Ve=vie? +ve? a er =] (77) 


This equation is not exact, but the author 
used it to calculate the approximate initial ve- 
locity of a star. 


6.3.2 Calculation of the Velocities of Stars 
near the Bursting Point from the Dots 
on Photographs 


The first dots in the photograph were not 
captured clearly by the revolving shutter. There- 
fore, the first dot of every trajectory was omit- 
ted. The difference between the coordinates of 
X4— Xz and that between ys — y2 were calculated 
and multiplied by the shutter constant of 12.5. 
From this we obtained the velocity of the third 
dot head x3,y3. Table 21 shows the data ob- 
tained by this method. 


Page 94 Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. 


Selected Pyrotechnic Publications of Dr. Takeo Shimizu, Part 4. Page 95 


Star 
No. (m) 
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Star 
No. (m) 
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Star Star 
No. (m) No. (m) 
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6.3.3 Calculation of Initial Velocities at the 
Bursting Point from the nearest Dots 
of the Trajectories 


The velocities v, and v, and the points x3 and 
y3 are obtained from Table 21. When these data 
are introduced into equation 77 using the value 
for k/p’, which is known, we can obtain the ap- 
proximate initial velocity of each star. 


When deriving equation 77, we assumed n = 
2 (see equations 74-1 and 74-2). Therefore, the 
value of k will not be the same as it was previ- 
ously when n = 1.93. Figure 26 shows log v, = 
1.218 when n= 2. 


logk =2-nlogv, =2—2.0x1.218 = 1.564 
From this value we obtain 
k =0.366 


Therefore, the values of the stars Bg, By and 
G,, are calculated as 


£800 6 0366 
p 10.0 

Bae =n 0nss 
p 142 

ee = 90083 
p 145 


The units for the above values of the sectional 
density p' are kg/m’. In this calculation the 
value o was an average value of the initial and 
those at the point (x3, y3). The mass m was as- 
sumet to be contant during the initial period 
from the bursting to the dotted position of each 
trajectory of a star. 


The results of the calculations are shown in 
Table 22, where the symbols are as follows: 
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V = initial velocity (m/s) for each star 

V = average initial velocity for one star 

V = ratio of the deviations of the initial 
velocity and the average velocity 
for one star 

r/V = ratio of the probability deviations 

of the initial velocities to the 


average for a star 


pee 
n-l 


where n means the number of stars calculated 
for the initial velocity for one star. 


@ = angle of elevation of the star 


Table 22. Velocities of Stars and the 
Deviations of Velocity from the Mean Value 
and the Angle of Elevation of Each Star. 


Av/V 
(%) 
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Table 23. The Average Initial Velocity V , the Ratio of the Probable Deviation of the Initial 
Velocity R/V and the Angle Deviations from Vertical Plane of the Trajectories of the Stars. 


1 
2] DiBeHWie | 84.14 | 66 | +6 | +14 


oN |e aBePMe o)  ed 


17 | DyBpsWe | 36.84 | 42 [ — | 
LB | DaBi2SWie | 43.32 [7.6 | 0 12 
vt 4BrgSMg | 33.82] 4.0 
Ot D4BigSMg | 65.06 | 2:6 | 10 | +0 
Le D4BizPWe | A787 | 42 | tO 6 

ee md 4Bi2PWie | 68.06 | ST. 


26] DsBeWie | 97.42 | 88 | 2 | sa 
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___DsBePMz | 122.49 | 10.7 | +0 | +0 
vedSBiZHWe | TAT | ee 
von sBizhWig | 93.98 | 6.6 8 | tO 
summa SBIZHMg | 67-39 | 9-5. 

vo sBigHMg | 105.32 | 4.9 | 


Ne cette aa tanec scincucce acs ccna taateal cattcnbacteente 


aed DsB,2S 


_DsBePWe_ | 63.46 | 42 | 20 | 30 


_DsBi2PWie | 


DsByPMa | 77.25 | 29 | +0 | +0 
DsBi2PMg | 98.58 J 7-6 | 


_f...O5C12G12B12PWe | 
vmmn5G12PWe | ven 2 P ede seem Do 
0%: 
DsC12PWs 


* In the case of star C the image on the photograph did not show any dots because of the weak light intensity 
and the production of trailing sparks. Therefore, with the C stars, the author could not use this method to ob- 


tain data for the velocity of these stars. 


6.4 Calculation of the Initial Velocity 
of Stars using a 16-mm Movie 
Camera and Comparison of 
those Results with Those 
Previously Calculated 


The images from the 16-mm movie camera 
appear more clearly at the bursting point than 
images captured by the revolving shutter, espe- 
cially for C stars. However, the problems with 
images from the 16-mm camera are: the scale 
was different throughout the picture and the 
number of frames per second was not constant. 
Therefore the 16-mm movie camera was not 
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useful to determine the exact initial velocity of 
stars; only the data from the revolving shutter 
were useful. However, the 16-mm camera gave 
clear images of the burning of the burst charge, 
which was very useful in studying the action of 
the bursting charge. (See section 6.5.3.) 
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6.4.1 The Scale of the Photograph Using the 
16-mm Movie Camera 


On the graph of trajectories (Sect. 6.2.1), the 
scale of the images was investigated. The ar- 
rangement of the star positions is shown in Fig- 
ure 31. The deviation of the flight surface of the 
photograph is ignored, and the surface of the 
flight paths of the stars is denoted by XYX'. 
The stars that flew from the bursting point flew 
along this surface. From the position of camera 
Z, the flight paths are photographed not on the 
real surface, but as they are on an imagined sur- 
face that would be perpendicular to the axis of 
the camera. The camera was set so that its opti- 
cal axis came perpendicular to the ground, and 
the flight surface of the stars also came to the 
ground in the same way (i.e., perpendicular to 
the ground). Therefore, the scale of the image 


Flight Surface 
on Photo 


Flight Surface 
of Star ———» 


Bursting Point 


along a vertical line does not change with 
changes in position. However, the scale along a 
horizontal line changes according to its position 
(Figure 32). 


Table 24. The Difference in Distance from 
the Bursting Point Due to the Scale Effect. 


(m) 


The quadrant of the image of the flight sur- 
face Oec' (Figure 31) is the object of analysis. 
The distance ratio X/X, at the point X, = 40 m 
(40 m signal light) was calculated as 1.04 (Ta- 
ble 24). This means that in the image of the 
photograph, the scale is almost uniform every- 
where in the quadrant. 


Signal Light 


aie) 
Camera 


Figure 31. Arrangement of testing devices when using the 16-mm camera. 
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The symbols in Figure 32 show: 
Z_ Position of camera 
X-X'-X' Flying star surface perpendicular to 
the ground 
A Point opposite the camera axis 
(center of photograph) 
AZ_ Light axis of camera 
X'-0-0'-X Horizontal base line, including burst- 
ing point to 20 m signal light, 40 m 
signal light (the camera is aimed at 
the 40 m signal light so that the 
camera axis o'-X crossed X-X' at the 
point o' of X-X' making a right angle 
with it. 
Ao Difference between the altitude of 
the camera and the bursting point 
a-b-e-c Direction of flight surface. The stars 
fly along this surface. 
a'-b'-e-c' Area of the image of stars photo- 
graphed 
o" The position of the 20 m standard 
light on the photograph 


6.4.2 Number of Frames for Each Star 
Image 


The number of frames per second of the 
movie camera was not constant for every star. 
Therefore, it is necessary to obtain the number 
star by star. For this purpose, the number of 
frames was determined based on the time pro- 
gression of the image from the revolving shut- 
ter. This was possible because each star was 
photographed by both the movie camera and the 
revolving shutter camera at the same time. 


For the standard, two or three stars were se- 
lected so that they appeared in the image, if 
possible, parallel to the vertical line BAB' (Fig- 
ure 31). 


The method used to obtain the relationship 
of the number of frames and time is as follows. 
First, from the data of the revolving shutter 
camera, the relationship of the flight distance 
from the bursting point and the time is obtained. 
Second, from the data of the movie camera, the 
relationship of the flight distance and number of 
frames is obtained. With the combination of the 


The Direction of 
the Camera Plate 


Flight Surface 
of Stars 


! 
:20m Xo 5 40m , 


/ 


Figure 32. The relation of the flight surface of stars with relationship to the ground and that of the 


image on the photograph. 
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two we have the relationship of the number of 
frames to the distance. 


The graph in Figure 33 shows the procedure 
for this method. From the image of the revolv- 
ing shutter camera, the star images (dots) with 
their relationship to time and distance are plot- 
ted in curve B against the distance scale (A). 
Then from the image of the movie camera, the 
closed points in the star light line are plotted as 
(C) against the distance scale (A). The relation- 
ship between (B) and (C) is obtained through 
the medium of (A). The number of frames per 
second was calculated as 80.6 as indicated in 
Figure 33. 


6.4.3 Comparison of the Data of Initial 
Velocities of Star from the Revolving 
Shutter Camera and the Movie 
Camera 


Using the above method, the initial velocity 
within the 4 circle of each image from the 
movie camera can be obtained. The process is 
as follows: The relationship of the number of 
the frame (7) and the flight distance of the star 
(D,») is plotted in a graph. The curve is extrapo- 
lated to the bursting point to obtain the relation- 
ship D,,:n. This then is converted to the initial 
velocity (dm/s). 


The data for 8 mm stars were obtained by 
averaging the flight paths. The data for 12 mm 
stars were obtained from each flight path then 
averaged. 


The 8 mm stars did not produce good im- 


0.1245 = 


—-_ 6.2 cuts x 0.02 s — 


(B) , 0,02 s 


(A) > 
30m | 
I 


ee ee eee 


ages because of the weak light from the prime 
composition, so the data were not very accurate. 
Table 25 shows the results. 


Table 25. Comparison of Initial Velocity 
Data of from the Revolving Shutter and the 
Movie Camera. 


Note: V' =data from 16-mm movie camera. 


Shell 28 showed a very large difference be- 
tween V’ (the movie camera) and V (the high 
speed camera). This may be due to the differ- 
ence in placement of the movie camera, which 
was set on the second day. The data from Shells 


From image of the flight path of 
a star cut by the revolving shutter 
from the bursting point. 


Distance from the 


aon I bursting point 


(ores 45 50 


<—— 10 frames ——— 


5 60 


Number of shutters 
from burst time by 
the movie camera 
from the burst point 


The frames per second = 10 frames/0.124 s = 80.6 


Figure 33. Calculation of the number of frames per second from sample shell star 491. 
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34, 44, 45, and 47 coincide very well for both 
cameras. The data for Shells 34, 42, and 43 are 
quite close for both cameras. 


6.5 Initial Velocity of Star Due To the 
Explosion (Bursting) 


Using the data for the initial velocity at the 
bursting (explosion) of the shell, a formula for 
calculating the initial velocity of the stars was 
developed based upon the type of shell con- 
struction. 


Generally, there were fairly large variations 
in the manufacturing processes, especially with 
different pasting operations. Therefore, when 
building the empirical formula, we must include 
theoretical considerations to avoid misleading 
results. 


6.5.1 The Influence of Various 
Manufacturing Processes on the 
Value of Initial Velocity of the Star 


Before this work, we must first examine 
how the initial velocities vary even when we 
use shells of the same construction. The sample 
shells were of different types. In this study, the 
author did not use shells of the same construc- 
tion; none the less, the analysis could still be 
made provided that appropriate theoretical con- 
siderations were taken into account. 


Shells that could be compared: 
|__| No. | Description 
49 


Ds5S12G12B12SWs 
D;G12SWs 


D;B12SWs 
DsS12G12B,2PWe, 
DsG12PWe 
D;5B12PWs 


The initial velocities of the above shells 
from Table 22 are as follows: 
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(i) For Shells 49 and 50, comparing stars G12 


Shell 49 Shell 50 


Star V Star V 
No. (m/s) No. (m/s) 
501 
502 
503 


504 
505 
506 
507 
508 
Average 


(11) For Shells 49 and 41, comparing stars Biz 


Shell 49 Shell 41 


Star V Star V 
No. (m/s) No. (m/s) 


Shell 53 Shell 54 


Star V Star V 
No. (m/s) No. (m/s) 


Type PEE ee 


DsBeSWe | OF 


(DsGiSWe bee it 102 


DiBrSWs | 59 oe ee oo mee re 225. 


D.Bi2SWs 
Note: W = the value of the bursting charge 


subscript o = full star shell 


(iv) For Shells 53 and 45, comparing stars B12 


Shell 53 Shell 45 


Star V Star V 
No. (m/s) No. (m/s) 


Average 


In the above four cases, the percentage of the 
average values are calculated as follows: 


51.8—56.0 


(1) ae 100=-8% for the stars Gi2 

(ii) ee 100=-8% for the stars Biz 

(iii) “ee 100=+1% for the stars Gy2 
65.0 — 63.5 


(iv) ————— x 100 =+2% for the stars By2 
65.0 


The above shows that the variation in the initial 
velocity due to various manufacturing processes 
can be as much as 8%. 
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6.5.2 Comparison between the Initial 
Velocity for Ring Shells and Full 
Star Shells 


Samples of full star shells were investigated 
using the movie camera. The flying stars were 
photographed as they were for the ring shells. 
Figures 34 and 35 show the construction of the 
image of a full star shell. (A) is the horizontal 
view, O is the bursting point of the shell, and 
O’z is the center of the axial line to the camera. 
In such an arrangement the star locus OP of 
distance x; gives the maximum flight distance 
Ox';, which appears on the image as line Ox. 
(B) is the vertical view from which we can ob- 
tain the maximum distance between the burst- 
ing point O and the highest point y;, which de- 
notes the image of a star OO. 


The values for maximum flight distance of 
stars obtained by the above method are shown 
in the table at the top of the previous page com- 
paring them with the ring star shells. 


To compare the flight distances of the full 
star and ring star, the formula: 


Ax x, —%, Ay = voi 
x x J J, 


was used to calculate the following: 


DsBgSWe | 0.353 | 0.512 


DsGy2SWe 
D,Bi2SWe 
DsBi2SWs 


The distance x,; was measured from the bursting 
point to the end of the blue (B) or green (G). 
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\ 
The maximum, 
vertical flight \ 
distance : 


The maximum \ 
horizontal flight ‘ 
distance \ 


Figure 35. Vertical view of the image from 
the movie camera. 


Figure 34. Horizontal view of the image from 
the movie camera. 


From the above results we know that gener- 6.5.3 Observation of the Burning of the 
ally the initial velocities of the full stars were Bursting Charge 
greater than those of the ring stars. This might 
be caused by the formation of a large resistance 
by the weight of the fully arranged stars in the 
case of the full star shell. Therefore, the effi- 
ciency of the bursting charge of the full star 
shell may be greater than that of the ring star 
shell. It must also be considered that the igni- 
tion powder (prime) on the stars other than the 
ring stars may also affect the initial velocity of 
the full stars. 


Figure 36 shows the results of some photo- 
graphs taken for the experiment through the re- 
volving shutter. The images show the part near 
the bursting point. 


For reference, the flight distance x,, and y,; 
of typical commercial shells are: 


Yor (m) 


From these results we know that much greater 
flight velocities might occur in typical shells 
than in the test shells used in the author’s work. 
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Figure 36. The images of the photograph 
through the revolving shutter. The images show 
the part near the bursting point. The photo- 
graphs are arranged in the following order: 


Ds5Bi2HWs Ds5Bi25Ws Ds5Bi2PWs 
DsBi2HWig = DsBr2SWie = 

DsBi2HM, DsBi2SM,4 DsBi2PM4 
D;5B,2HMs Ds5B,25M, Ds5B,2PM, 


In each photograph we observe a flash of 
light from the bursting charge, which was still 
burning after the shell broke. This burning is 
ineffective in propelling the stars because it 
cannot create pressure in the open air. It could 
be said that the larger the after-light, the smaller 
the efficiency of the bursting charge. (However, 
it may be necessary to investigate the light in- 
tensity of the bursting charge.) 


From this standpoint, using the same method 
of paper pasting, the efficiency of the bursting 
charge is the best with H, medium with S, and 
least with P. The P burst powder is the least 
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effective, and it was seen to burn in the air for 
about 7 meters in the case of D;B;2PWs3. How- 
ever, when the thickness of the paper pasting 
was greater, as in Ds5B;2PMsg, the flash of light 
was small, which shows that the efficiency of 
the bursting charge was greater. 


For comparison of the burning velocity after 
the shell break, movie camera photographs (A) 
Ds5B,2HWs, (B) DsB12SWs and (C) DsB12PWs 
are shown in Figure 37 (A—-C). Each of the pho- 
tographs were arranged on every two frames. 
The time interval for the two frames is 1/60 
second. 


Figure 37 (A). Ds5By2HWs with movie camera. 
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D5B142PWs with movie camera. 


Figure 37 (C) 
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6.5.4 The Equations for Calculating the 
Approximate Initial Velocity of a Star 


This equation is for one of many stars in a 
chrysanthemum type shell. Then we have from 
equation 31: 


-1/2 


V=K-E-R' (m/c) (78) 


The value of E depends upon fw from equa- 
tion 30, and the value of R depends on the 
number of pasting layers on the shell, the qual- 
ity of the pasted paper, and the burn velocity or 
the vivacity of the bursting charge. The stan- 
dard paper strength was 230 kg/cm’ for the 
pasted warp and woof and the standard thick- 
ness is 0.08 mm per sheet. 


The following symbols are used: 

nthe number of layers of the standard Japa- 

nese paper pasted on the shell per sun 

the number of layers of other type paper 
pasted on the shell per sun 

a comparison of the strength of number of 
layers n of the standard paper with that of 
the other paper n' when the pasted layers 
have the same strength (i.e., a = n/n’) 

A the vivacity of the bursting powder, which 
is denoted by the reciprocal of the total 
burning time T in the air (i.e., n = 1/7, 
where T is in seconds) 

a exponent for the sectional density of a star 

2 exponent for the total energy of the 
bursting powder 

kn aconstant concerning the efficiency of 
acceleration on the velocity of a star 

N  N=knl/A and could be called “specific 
pasting number” 


Then with £& and R we tried to assume that: 


E" =e( fir) =( fay’ (79) 


kn N 


RP =R(n,a,A)=e “=e @ (80) 


For equation 80, the following conditions were 
fulfilled: 
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an > © R= $1 
an > 0 R= 0 
A> o R= 1 
A > 0 R= 0 


Only two or three data were obtained from 
experiments; however, when conditions for 0 or 
co are included, the number of data increases to 
four or five. Including the latter data increases 
the accuracy of the experiment. When equa- 
tions 79 or 80 are substituted into equation 78, 
and the value of the exponent 4 for m/o is 
changed to a more general value a, we have 


V = K(f) (m/o) em (81) 


Then the unknown values a, B, a, N, and K are 
determined from the experimental results. 


(1) Determination of a 


The value of o is determined from the ratio of 
the initial velocity of stars of two types of shells 
that have the same manufacturing parameters 
except for the size of their stars. Namely with 
the size of stars Bg and Byp2 the ratio is intro- 
duced as: 


K(m)(™) ew [m | 


Ve. oO /p, on, 
v, aN | m (82) 
By K(fiv) m em —- 
Oo B,, O Bp 


From the logarithm of equation 81, we obtain: 


V. 
log— 
a= (82') 


. 


Introducing a symbol V for the average value 
V of the initial velocity in Table 23, we have 
the following relations: 
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D,BgHW 16 
D4Bi2HW1¢ 
D,BgHM, 
D,B,2HM, 
D,BgHMg, 
D,B,2HMs 
D,BgSWs 
D,Bi2SWe 
D,BgSWi, 
D4Bi2SWi6 
D,BgSM, 
D,Bi2SM4 
D,BsSMz, 
D,Bi2SMs 
D,BgPWs 
D,By2PWe 
D,BgPW, 
D4B12PWi6 
D,BgPM, 
D,By2PM, 
DsBgHW 16 
Ds5B12HW1¢ 
D;5BgHM, 
DsBi2HM, 


D;BgSWs 
D;Bi2SWs 
D;BgSWi6 
DsB12SWi6 


D;BsPWs, 
DsB12PWs 
DsBgPM, 
DsBi2PM4 
D;BgPMs 
DsB12PM, 


The * indicates an irregular value. After re- 
jecting the values with an *, the values for the 
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remaining 16 data points were averaged. The 
values for the sectional density were taken from 
p, (Section 4.3.1 —(5)). 


WAR _ Po, _ 0.92 _ 
(™/ ),. Py, (1.28 


Therefore, from equation 82: 


es log 1.045 
~ log0.719 


=-0.133 


(2) Determination of N (Using Japanese 
Paper as the Standard) 


Initial velocities were compared for shells 
with the same construction but with a different 
number of pasting layers using the standard 
Japanese paper. From equation 81, the ratio 
Vy, and Vy, 1s denoted as: 


16 V. 
N=- xlog} “6 83 
0.4343 e( “fy | ? 


For each type of bursting charge the values of 


yf were calculated: 
Vw, 


Burst Charge H: 


D,B,2HW 4, 
D,Bi2HWe 
D;5BgHW,¢ 
D;BgHWe. 


Average value for (H) 1.179 
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Burst Charge S: 


D,BgSW 6 
D,BgSWs 


2 
D;BgSWi6 
D;BgSWs 
DsB12SW 16 
DsBi2SWe 


D,BgPW 6 

D,BsPWs 

D4B,2PW 16 

D,B.2PW, 

DsBeP Wie 1.801 
D;BsPWs 


Average value for (P) 1.519 


Therefore the value of N for each bursting 
charge is obtained from equation 83: 


For potassium chlorate bursting charge (H): 


N= 2 x log1.179 = 2.633 
0.4343 


For Black Powder bursting charge (S): 


N, = - x log1.223 = 3.218 
0.4343 
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For potassium perchlorate bursting charge (P): 


Np= = x log 1.519 = 6.687 
0.4343 


(3) Determination of a 


The value for a was determined by compar- 
ing the velocity of the stars that were from 
shells that have the same construction except 
for the type of pasted paper but with the same 
total strength: 


For the relationship W4g—Ms: 


17 1 1 
Wiis Sea oe) 
i 
Therefore: 
V. 
lop —* -oasasy{ 1) (84) 
Pric 8a 16 


Similarly for We—M, 


=0.4343N e - ‘| (85) 
M, 4a 16 
From equation 84: 
ne Se (36) 
Ba yp 
. Vay, : 1 
0.4343N 16 
Similarly from equation 85: 
Ge , = : (87) 
4 Vi 
log} "> 
Via 4“ 1 
0.4343N 8 
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From equations 86 and 87 we have the following values of a: 


oe 


D;B,2PM, 
Average | 2.331 


From the above calculations we have a = too large. Therefore, the values were averaged 
2.33 as the average value. The values marked from the remaining 14 values of a. 
with * were rejected because the deviation was 
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(4) Determination of the value B 


B was calculated by comparing the initial 
velocities of 4 and 5 sun shells. 


From equation 81 we have: 


Taking its logarithm, we have: 


poke |e ay | 


By using the same type of bursting charge, fis 
eliminated as follows: 


Vp, Wp 
rule] 


For each type of bursting charge the follow- 
ing values were calculated: 


Ww 

H: —= =0.234/0.107 = 2.19 
Wp, 
wy 

S: + = 0.232 /0.106 =2.19 
Wp, 
Ww 

P: —* =0.256/0.116 =2.21 
Wp, 


7 
And the values of > were calculated as: 
Dg 


Burst Charge H: 
DsBsHWii¢ 


D,BsHWi6 


DsBi2HW46 
D4Bi2HW46 


DsBgSW16 
DsBsSWi6 


DsB125Wes 
DsBi2SWe 
oe 


1.521 
0.8.8 : 


D.B..SM, 
D4B125M4 
DsBi2SMgz 
D,B,2SMg, 


1.812 


Shell V 
Symbol 
(m/s) 


DsBsPW 16 
D,BsPW16 


Average value for ( 1.513 
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The values with * are rejected because the 
deviation was too large. The value of B for each 
type of bursting charge is as follows: 


H: £=log1.325/log2.19=0.357 

S: B=log1.446/log2.19=0.470 

P: B=log!.513/log2.21=0.552 
Average 0.450 


(5) Determination of the Vivacity A for Each 
Bursting Charge 


From the definition in 6.5.4 
N= a therefore A= uals (89) 
A N 


Therefore, when the value of kn is known, we 
can obtain the value of A. Alternatively, the 
values of A are obtained by taking the recipro- 
cal of the burning time in the air from Table 13. 


A, =1/0.167 =5.99 (1/s) 
A, =1/0.240 =4.17 (1/s) (90) 
A, =1/0.518 =1.93 (1/s) 


Editor’s note: The symbol “As” used here for the 
vivacity of Black Powder bursting charge must not 
be confused with the same symbol used elsewhere in 
connection with the ballistics of flying stars. 


When the value of As is introduced into 
equation 89, the value of kn is obtained as: 


kn = A,N, =4.17x3.218 =13.42 


This value for kn and the values for Nj and Np 
are substituted into equation 89 to obtain the 
values for 4;; and Ap relative to the vivacity of 
As: 


A, =kn/ Ny, =13.42/2.633 =5.10 


(91) 
A, = kn! N, =13.42/6.687 =2.01 


When comparing equations90 and 91, we 
know that both values nearly coincide. This 
shows that the values of A from the movie cam- 
era taken of the burning bursting powder on the 
ground and those taken from the explosion of 
the bursting charge of the shell coincide fairly 
well. (However in this case, the value of As is 
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used commonly for both cases.) From here on, 
the value of A from the explosion of the shell is 
used. 


6.5.5 Adaptability of the Empirical 
Equation of Initial Velocity of a Star 


The coefficients in equation 81 were all ob- 
tained as described above, and hereafter the 
author examines how they are adaptable to the 
empirical initial velocities again. 


(1) Calculation of K 
From equation 81, as the first estimate, the 


following calculation was performed: 


K'=V (fw) ” (2), Len (92) 


Oo 


For the values of f, the ones calculated using 
equation 4 from 4.3.2 were used: 


fy =9.736x10° dm =7.36x10* m 

fp =0.712x10° dm =7.12x10* m 

fg =0.287x10° dm = 2.87x10* m 
To know if these values for fare reasonable or 
not, we have to see whether or not the calcu- 
lated values of K are independent of the values 
of f; For the values of m/o, the values of p; 
from Table 5 were used: 


(=) =9.2 kg/m? 
Oo)». 


(=) =12.8 kg/m? 
o By 


The other constants came from 6.5.4: 


a = 0.133 
B = -0.450 
Nu = 2.633 
Ns = 3.218 
Np= 6.687 


a= 1.00 (standard Japanese paper) 
a= 2.33 (kraft paper) 
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The results of the calculation are as follows: 


Burst Shell 
Powder : log K’ 


The symbol * shown in 6.5.4 means rejected once. 


The symbol ** means rejected twice. For the calculation of K’, those with ** were rejected to minimize the er- 
rors. 


(15)* was also rejected because the result was not reasonable because it produced the relationship V(15)<V(i4), 
which is a contradiction to general physical law. The average values of log K' for the types of burst powder are 
shown in Figure 38. From Figure 38 we see that the values of log K’ are grouped into four areas: BgDy, By2D4, 
BsDs, Bi2Ds. For each group, the values of K decrease in the order of P, S, H. The position of each center 
somewhat deviates. The sectional density may cause this. And the difference between Kp, Ks, and Ky, would 
come from the difference in the value of vivacity rather than the value off Therefore, first the position of the 
center of log K is corrected. Next, the difference of log Kp', log Ks', and log K;', which are the result of the influ- 
ence of vivacity, are rejected values, and finally, the constants are determined. 


log Ky, log K;, 


H 0.2668 0.1809 
D,38 0.3727 0.2856 
P 0.4151 0.3955 
H 0.2413 0.2384 
D,3S 0.3139 0.3284 
P 0.4306 0.4015 
Average 0.3401 0.3051 
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Figure 38. Average values of log K"'. 


Using the derived constant K", we have: 
(93) 


In these equations, Ke and ae are the av- 


erage values of K for each case, and a' is a de- 
rived exponent. Both terms in equation 93 are 
divided by each other and taking the logarithm 
we have: 

log Ky, - log K,,, 


, 


© toe), /lm), | 


(94) 


0.7 


0.6 


log K" 


0.5 


0.4 


in place of log Ki and log Ky , using the av- 
erage value of log K;, and log K;,, and intro- 
ducing the values of YA we have 


1 0.3401-—0.3051 __0.244 
log(9.2/12.8) 


And from equation 93 we have: 
log K" = log Ky, = log). 


= log Ke - log (™/) 
In this formula 


log(™/ ) = ~0.244 log 9.2 = -0.2352 
log(™/ ): = ~0.244 log 12.8 = —0.2702 


These values were introduced into the above 
equations: 


log K" =log Kj, + 0.2352 =log Kk, + 0.2702 
(95) 


For each value of the log K’, they are corrected 
with equation 96. 


log K" =log Ky —0.2352= log Ky, — 0.2702 
(96) 


The next table shows the results from equa- 
tion 96. The average values of log K” are 
shown in Figure 39. 


Group Interval 


Figure 39. Average values of log K". 
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From Figure 39 we see that the influences of 
sectional densities are mostly rejected. 


Second, from the values K” the influence of 
the vivacities were rejected. When K" is de- 
noted as: 


K"=KA* (96) 
Then for each bursting charge we have: 

KU =KA,, KH KA, Ke KA, (97) 
Taking the ratios: 


AUR. oie 
And taking the logarithms of above: 
log K;, —log K; = x(log A,, — log A,) 
log K{ —log K; = x(log A, — log A,) 
From above: 
_ log K;, — log K; 
log A,, — log Ap 
_ log Ks; —log kK; 
log A, —log A, 


(98) 
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The average values of K" for H, S, and P are obtained as follows: 


Average 0.4836 || Average 


From equation 90 the following are obtained: 
log A,, = 1og5.10 = 0.7076 
log A, = log 4.17 = 0.6201 
log A, = log2.01=0.3032 


Introducing the above values into equa- 
tion 98, the value of x is obtained. 


ve 0.4836 = 0.6637 _ 
0.7076 — 0.3032 
ve 05815 = 0.6637 _ 


~ 0.6201—0.3032 _ 
Average —0.352 


0.445 


0.259 


From equation 97: 

log K = log K;, —xlog A, 
= log Kf —xlog A, 
= K, —xlog A, 


The value of K, which is amended by equa- 
tion 97, is thought to be the value of K that is 
independent of the vivacities of the bursting 
charges. The values of K amended by the values 
of K" from the above equations are shown in 
Table 26, with the following calculations: 
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Average 


xlog A,, =—0.352 x 0.7076 = —0.2491 

xlog A, =—0.352 x 0.6201 = —0.2183 

xlog A, =—0.352 x 0.3032 =—0.1067 
The average of the values in Table 26 is 5.91. 
The probability deviation of K is calculated as: 


—\2 
=(K -K) 
Y, =0.6745,{ 
He 
= 0.6745, ans = 0.496 
37 
Therefore, 


K =5.91£0.496 =5.91(1+ 0.084) 


That is, the value of K has an 8% probability 
deviation. As a test, we examined the average 
values of K for each bursting charge as follows: 
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Table 26. Proportional Constants in the Empirical Equation for Initial Velocity of Each Shell. 


a Shell Shell Shell 
No. No. No. 


1 
a 
Fic 
8 


kn 


These data are illustrated in Figure 40. _ in 
anA (99) 


In Figure 40, the difference by fis not clear. 
Therefore the value of f calculated using equa- 


ian oF 43.2 ie-cottect and weetal. Introducing the following values that were 


already obtained: 
(2) Determination of the Empirical Equation 


for the Initial Velocity of the Star a =—0.133, a! =—0.244, a +a" =—0.377 
B=-0.450, x =—0.352, kn =13.42 


K =5.91(1+0.084) 


The tentative equation 81 is rearranged as: 


a kn 
me 2 
=K(f a)” (“) em we have the following equation for practical use: 


—0.377 13.42 
The final equation for the initial velocity of V =5.91(1+0.084)( fo) (“) Ae @ and 
the star is denoted as follows, adding the condi- Oo 

tions of equations 93 and 97: (100) 


0.9 


0.8 


log K 


0.7 


0.6 


Group Intervals 


Figure 40. The constants K that are independent of K. 
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MKS units should be used in equation 100. Equation 31 | Equation 100 


When comparing equation 100 with equa- 
tion 78, which approximates the equation for 
initial star velocity, we notice the following: 


The exponents roughly coincide with each other. 
This shows that the idea of the development of 
the empirical equation of the initial velocity did 
not deviate very much. 
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Chapter 7 — Equations for Practical Use in Designing 


Shells 


7.1 The Characteristics of Stars and the 
Preparation of Tables of Star Ballistics 


Equations 71 and 72 are written again: 


-1.075 


A l-z') +b 
ron [4 ioe! rvs (71) 


(72) 
When ballistic functions are set as follows: 


F. =F. (A,,7’,b) 


Lhad: Gey abl “Ot 
-{ [4 og ea 


F. = FASEB) 


(I \ b —1.075 
_ ra Ay —E i 
=| f [ 3 }toe is dt 


(102) 

equations 71 and 72 become: 
y= Vol", (103) 
aay iw oe (104) 


The tables in the Appendix list the calcu- 
lated values from equations 101 and 102 where 
b=0Oandb=0.5. 


7.2 Investigation of Star Ballistics in the Air 
from the Ballistic Tables 


From equations 103 and 104 we have: 
v Pe: 
see Sania 


Therefore, for a single colored star, the function 
F, shows the ratio of the flight velocity to the 
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initial velocity, and F;. shows the ratio of the 
actual flight distance to the imaginary flight 
distance with no air resistance. The larger the 
values of F,, and F’, are, the greater the move- 
ment of the star. Figures 41A and B and 42A 
and B show the values of the functions F, and 
F,, when b = 0.0 and b= 0.5. 


The values of A, are denoted as follows from 
equation 70' without the subscript 1 for R, 7, 
and 6 

Uf 
R/T’ 6 


A, = 0.774)" 


In the above equation, R/T’ means the burning 
velocity and 6 means the density of the burning 
layer. Therefore, when the value of v, is con- 
stant, the larger the burning velocity, the 
smaller the value of A,, and the larger the spe- 
cific density, the smaller the value of A,. Fig- 
ures 41 and 42 show that the smaller the value 
of A,, the larger the values of F, and F,.. There- 
fore, when the burning velocity R/T’ and the 
specific density of the star are larger, then the 
star flies further. 


The value of b is denoted from equation 67' 
as follows: 


3 
| Ba) | Se Oa ly 
R, 6, 6, 
3 3 
R, 6, 6, fs R, 6, 1 
R 6, 6, R 6, 


If there is a relationship between the layers of a 


star: 


O&O S05 Sats 


the value of b becomes larger and larger as 
burning proceeds. Comparing b = 0 and b = 0.5 
in the graphs in Figures 41 and 42, we see that 
the decrease in velocity of the star is less with b 
= 0.5 than with b = 0.0. Furthermore, the flight 
distance for b = 0 is greater than for b = 0.5. In 
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Figure 41. Velocity functions. 
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a 


Figure 42. Distance functions. 


other words, a star, of which the density becomes 
greater and greater, produces the best effect. 
Special attention should be paid to the condi- 
tion when 5 = 0 and F’, = 0. However, when b > 
0 under t= 1, F, has a limited value. 


7.3 The Method For Obtaining the Function 
F,(7, A,, b) from F,(7’, A,, 0) 


The ballistic tables in the Appendix were 
made for b = 0 and b = 0.5. For other values of 
b, they do not correspond. However, with the 
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25 
30 


0.0 
00 01 02 0304 05 06 07 08 09 1.0 


a 


0.0 01 02 0304 05 06 07 08 
q' 


following method, it is possible to obtain the 
data from F,(t', A,, 0). The function of the prob- 
lem is set as: 


PAK A,, b) 
It is equated with: 
F,(¢,,4,,0) = F,(¢', 4,5) (105) 


Both sides are substituted into equation 72 to 
give: 
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A oe —1.075 
{1-(4)ios(1-«:) 
(4); (27) ) 
= —_|—_— og ——_ 
cS ace ee 


From this equation we have: 


1x9 
log(1-r,) = log) +8 

a u (106) 
7, =1=3 eau? a =t(r',b) 


Therefore, the value of F,(7,A;,b) is equal to the 
value F',(7,,45,0), which is calculated with the 
value 7. Therefore from the table for F\(7,A,,0) 
we can obtain the values of F\(7,A,,b). For this 
purpose a table of F(7,b) was prepared. The 
values of F,(7,4;,b) are simply not available. In 
this case it is necessary to calculate the function 
F(7,A;,b) and to integrate step by step. 


7.4 An Examination of the Equations of 
Initial Velocity of Stars 


From equation 100, setting three equations 


as: 
7 = (fo) 
0.377 
f= (2) (107) 
(oy 
f = pice 
we have 
V =5.91x(140.084) fA -f (108) 


Here it should be noted that equations 100 and 
108 are empirical and apply for loading densi- 
ties in the range 0.27 ~ 0.33. 


Special attention should be given to Figure 45 
with the function f;. When different types of 
paper are used with the same value of an in to- 
tal, the value of f; moves along the same an 
curve. When the value an is relatively small, 
the curve f; does not show the maximum. How- 
ever, when an shows a high value, the curve 
shows a maximum as an = 12 or 16. For exam- 
ple, when pasting with standard Japanese paper 
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250 


200 


150 


(fy) 


100 


50 


2 4 6 8 10 12 14 16 18 20 
f@ x 10-4 kg/m? 


Figure 43. Function f; for the initial velocity 
of star. 


on a Shell where an = 8 and A = 5.0, we have f; 
= 0.41. If the bursting charge is replaced with 
that which shows a slower burning rate (4 = 
1.5), the same value for f; as above can be ob- 
tained when an = 12. In this case, one should 
paste 12 layers per sun. When the pasting layers 
of paper increase to a large amount, the effect 
of the force on projecting stars by the bursting 


1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 


(f,) 


0.3 
0.2 
0.1 


0.0 
2 4 6 8 10 12 14 16 18 20 


p/o kg/m? 


Figure 44. Function f> for initial velocity of 
Star. 
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Figure 45. Function f; for initial velocity of star. 


charge increases in a similar manner. 


7.5 The Relationship between the 


A 
Deviations of Initial Velocity ae and 


Ax, 
the Deviations of Flight Distances —+ 
x 


The flight distance of a star along the hori- 
zontal axis during the burning of the first layer 
of the star is shown in equation 104. 


x a VIF. (4,,7',5, ) 
If the construction of the star is unchanged and 


T',, 7’, b; are constant, the equation is differen- 
tiated as: 


dx, dV dF. 


Bh es 109 
x V se 
On the other hand: 
dF. a A, OF. \ dA, (110) 
F. F, OA, A, 


If, as an approximation, the exponent of V is 
put as 1, 
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A, 0.7747 (111) 


WO; 


If this is then differentiated by the logarithmic 
method, we have: 

dA, _ dV 

A V 


Ss 


Introduce this into equation 110: 


dF. {| A, _OF, dV 
Ef ey 0A, V 


x 


(112) 


Introducing equation 112 into equation 109: 


PF 
dx, _ dV 4 be Fi (113) 
Rice oe F. OA, 
In equation 113 we set: 
A, OF. 
pd ate Se 54 (4,2',3) 114) 
F. OA, ial 


and when dx, and dV are replaced by the devia- 
tions Ax and AV, equation 113 becomes: 


calle 5 a (115) 
x V 


Equation 115 is generalized as: 
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Ax. Avy, 


Xx; Vo. 


(115') 


Equations 115 and 115' are equations to denote 
the relationship of the deviation of the initial 
velocity and horizontal flight distance. The val- 


ues of XA, =| 1+ BOE. are shown in the Ap- 
F OA 


Ss 


pendix in Table 2. 
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Generally, as the value of A, increases, the 
values of X A, become smaller. This shows that 
as the value of A, increases, small deviations in 


the value of the initial velocity do not have 
much affect upon the flight distances. 
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Chapter 8 — Examples of Designing Chrysanthemum 
Shells 


8.1 Comparison of Calculated and 
Experimental Values 


can be examined by the probability deviations 
of the constant n' and log v,. Therefore, the au- 
thor only compared the results obtained from 
the empirical equation with the results of ex- 
periment. 


The accuracy of the empirical equation for 
the initial velocity of a star was examined in 
Section 6.5.5. The accuracy of flight velocity 


Example: From the Data Given Below, Calculations of the Flight Distance of Stars. 


Construction D;G,.PW. 
ee 2) ee es <1, ae 


|...9.00055kg | 


nnBursting charge | | 


ae 


[1620 kg/m® | 
[1650 kg/m? | 
odd 004 SS ne 

wn;00250 mr. 


P 


o ee ee eae 


H 


1620 kg/m 
1650 kg/m | 
9.00600 afm 
7 a 
La: a 


a a 


1640 kg/m” | 
1650 kg/m? 
0:00600 m_ 
00250 ae 
9.00156 kg. 


P 


pe ASIN A tna NL | ome papel) . sapeese Pape | 
0.30 0.30 0.30 


Note: S was omitted because it is the weakest of the three and not useful except in some special cases. 


142x108 m_ 


(a) Calculation of Initial Velocity VK fh hs 


Equations: wien 
K=5.91 


f = (fo) 


o-(2 
oO 


_ 13.42 
f. = 4°36 anA 
3. 


(Table 7 in the Appendix) 
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Calculation Number 


3.14x(0.00600)" 3.14x(0.00610)? 
0.00147 0.00156 
13.00 13.35 


f; (Table 7 in Appendix) 


Multiplying K by the values of f), 2, and f; 
we have: 


Calculated V 


Experimental V 
Calc. VIExp. V 


The above table shows the coincidence of 
the values of initial velocity between the calcu- 
lated and experimental. 


(b) Calculation of the Flight Distance of the 
Star in the Horizontal Direction 


During the time from the start to the end of 
the burning of the blue or green zone, excluding 
the bright flash from the burst charge, the flight 
distance was calculated. For this process, the 
final velocity v; was first calculated. 


Equations. 
Were, (4,5,.7;) 
A, =0.774V°”° ail 
Ro, 


T. 
T/=—, 
i 
R, 
a: 
Z== 
T, 
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ee ae ae ee ae 
le OE a Eo OS 
0.00435 | 0.00600 | 0.00610 

1620 1620 1640 
41.97 53.99 45.96 
7.80 11.79 17.82 

OOO 2 eS |. 
0.425 0.583 0.583 


Assuming b = 0 and using the values of 4, and 


T,'1n the table of F., we have: 


Multiplying the above three values mutually, 
we have: 


The experimental values were: 


Experimental x; | 28.0m | 52.1m | 63.1m 


Star Number 


The value of 7; was measured experimentally. 
The result shows that the calculated values of 
the flight distance were 3.6 to 4.6 m smaller 
than the experimental values. 
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8.2 An Example of Designing a 
Chrysanthemum Shell 


In this paper only an example of designing 
calculations was introduced. In practical cases 
more important factors need to be considered: 
for example, ease of manufacture, efficiency of 
explosive, etc. 


Example: A Flower of 80 m Horizontal 
Radius, Amber to Blue Color with a 
Horizontal Flight Distance of 40 m Blue and 
40 m Amber. 


(a) Calculation of the Second Burning Period 
(Blue Colored). 


From the given condition: 
x, =40.0=1,T,F (A, .b,.7,) (i) 


In this equation, v, is the final velocity of 
the star at the end of the first period and at the 
same time it is the initial velocity of the second 
period. x2 is the distance traveled in the second 
period. The symbols for the first period are 
noted by a subscript of | and those of the sec- 
ond period are noted with a subscript 2. There- 
fore: 


T;=T, 
A, =0.774y,” a = 0.774." 7 
From experiment (Tables 5 and 16): 
w, =0.00195 m/s 
5, =1620 kg/m* 
Therefore 
A, =0.2450v,°° (ii) 


On the other hand the value of vo 95 is calculated 
as: 


Vous = VF, (A,, 00.95) (iii) 
From (ii) several v; values are obtained, and 


from (iii) v.95 values are obtained as shown in 
the following Table. 
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The relationship between the data is illustrated 
in Figure 46. 


Figure 46 shows that when the values of 
A, exceed 9, the final velocity of the star does 


not increase as much. Therefore, a value of less 
than 9 may be desirable. The larger the value of 
Vo.95, the less the star is affected by gravity. For 
this requirement, the value of A, should be as 


large as possible under the condition A, < 9. If 


we use A, = 9, then we have: 


v, =48.14, Vyo5 =3.13 


Then the radius of curvature of the trajectory of 
the star from equation 72' is: 
vy (3.13) | 


cos@ = = =-1.0m 
” z 9.280 


The value of 1.0 m is too small. This means that 
the final velocity of the star is too small. How- 
ever, if we possibly have stars with a large den- 
sity or large final velocity, the problem may be 
successfully solved. From (i) we have: 


ae. 
VF. (9,0,1) 


T= 


From the above Table we have: 
F. (9,0,1) =0.318 


Therefore, 
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4.0 


Ag 
44 12 13 14 15 16 
ate 9 po -e- 9 6-9-6 
3.0 5 9-0-8" e 
o7® 


Vo.95 (m/sec) 


50 


v4 (m/sec) 


Figure 46. The values that are calculated from v; and the function F,. 


100 


bee ee VF,(A, ,0,z/) =v, =48.14 (iv) 
48.14x0.318 
From this value, we can determine the radius of Ve, (4, ,0,7/) =x, =40 (v) 
the second layer at the beginning of the blue: hn 
R, = wl, - 1 
= 0.00195 x 2.61 A, =0.774V~ bea 
“1 
= 0.00509 m 1 
- =0.774°? —______ 
eeprene 0.00234 1480 
= 0.93 
(b) Calculation of the First Burning Period = 0.2235V 
(Amber Color) From this equation the relationship of A, and 
From an experiment: V is obtained as follows: 
6, = 1480 kg/m* 
— ae 47? fT 0.60 
= 0.00234 m/s {4.3.1 (4) or (6)} In this case the final value of 7’ is unknown. It 


Connecting this period with the second period: 


is assumed with a proper value of 7’ that the 
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0.20 
Figure 47. The relation of v; and T, . 


v, =VF,( A, ,0,7; 
ae Mal @ [=| =| sa 
_0 iy F 10 


62.45 


49.21 | 50.14 | 50.82 

44.99 | 45.72 | 46.36 | 47.01 | 47.58 
41.35 | 41.91 | 42.38 | 43.02 | 43.42 
38.15 | 38.56 | 39.05 | 37.53 | 39.96 


The above data are plotted as a graph in Fig- 
ure 47, and the final velocity of v; = 48.14 
m/sec is recorded as a dashed horizontal line, 
and the points of intersection give the value of 
z';, these denote the relations of values A, and 


Cas 
T° 
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0.13 014 0.15 


And then, 
tT, =1-R,/R, 


From this relation: 


R, 0.00509 
R, = ' = ' 
l-z, 1-7, 
On the other hand: 
' T, = R, R 


= 427.4R, 


1 


~T-R,/R, w, 0.00234 


Combining the above two equations: 


T'= one 
1-7, 
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0.16 0.17 0.18 0.19 0.20 


lows: 


Using these values and the value of V, the travel 
distance x, during the first period is calculated 
as follows: 


x; = VLE, (4,,0,7;) 
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These values are plotted in Figure 48. 


In Figure 48, the condition was given that x 
= 40 m. It is shown in Figure 48 as a horizon- 
tally dashed line, and when we read the value of 
t; at the points of intersection with the A, 


curves, we have: 


From the data of (I) and (II) we have the two 
curves on the graph in Figure 49. 


The point of intersection denotes the values of 
A, and V that satisfy the conditions given be- 


fore: 
A, = 26.12 
v', =0.175 


For the value of A, we have: 


V =167.0 m/s 


Therefore, we have: 


28 


27 


0) 
24 


R= ue : 
l-f, 
_ 0.00509 
~ 1-0.175 
= 0.00617 m 


=6.17 mm 


Therefore the diameter of the star is: 


d =0.00617 x2 
= 0.01234 m 


=12.3 mm 


Then we have: 
LR 
= Fe 
_ 0.00617 — 0.00509 


0.00234 
=0.46s 


ihe 


The solution to the first period is completed 
with the above calculations. 


The burning time of the star in total: 
T=7,+T, 
=0.46+ 2.61 
=3.07s 


(II) 


0.18 0.19 0.20 0.21 


Figure 49. Graph for deciding the value of A, and t',. 
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The mass of the single star: 
m= Sa {(R —R})6, + R36;} 


0.00617) —(0.00509)° 
Se ( ) ( | 


x 1480 + (0.00509) x 1620 
= 0.00153 kg 
=1.53 g 


0.00153 
3.14x (0.00617) 
=12.80 kg/m? 


m 
Oo 


In the above calculation it was assumed b; = 0. 
This is checked: 


3 
(oa 
1 1 
i 1 boast ) 
1480 /\ 0.00617 
= 0.053 
This value does not disturb the condition b = 0. 


The above calculation concerns the design of 
a star flying horizontally, but the falling height 
of the trajectory is as follows, ignoring air resis- 
tance: 


1 
=— eT” 
Ra 78 


=5*9.80x(3.07) 
= 46.2 m 


The Design of the Shell 


First an assumption is made that there is no 
difference between ring star and full star shells 
with regard to the design conditions. As de- 
scribed earlier, the initial velocity of the full 
star is somewhat greater than that of the ring 
star. However, this does not affect our manufac- 
turing process. 


The initial velocity from equation 108 is: 


VSR “ff. f, 
=167.0 m/s 
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From this relation we have the value of /;: 
(vii) 


On the other hand: 


nly” 


2 (12.80) 
= 0.382 


When we use the bursting charge P, the charac- 
teristic values are: 


f =7.12x10*, 4=1.93 
For shell pasting, when we use Kraft paper 
(cement bag) for 6 layers: 
an = 2.33x6 
= 13.98 


From the values of an and A, the values of f; 
are calculated (or from Table 7, value f3) 
13.42 


_ 470.352 and 
re -A em 


13.42 
= (1.93) °° e@ 13.981. 
= 0.483 
Therefore, from (vii) 
167.0 
5.91x 0.382 x 0.483 
=153.1 


= (fa ‘ae 


From above, the value of f@ is 


i 


fao=71710 kg-m 


When we divide fa by f= 71200 m, we have 
1.007 kg. 


The quantity of cottonseeds a’ , which fixes 
the bursting charge @ , assumes w/a’ = 16: 


ow = w/1.600 
= 0.629 kg 
a’ = 0.629 kg 
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Therefore the total mass of the burst charge is: 
W=a0+a' 
= 1.007 + 0.629 
= 1.636 
=1.600 kg 


The apparent specific gravity of the bursting 
charge is 0.556 and the volume of the charge is: 


_ 1,636 


y= 
0.556 
= 2.942 liters 


When the volume is spherical, the radius is: 


1 
Re (2. tes 2942) 
4a 


= 0.889 dm 
=89 mm 
Including the radius of the star: 
R+d=101mm 
Therefore the inside diameter of the shell is d; = 
101 x 2 = 202 mm. For the 7 sun (~8 inch) 


shell, the recommended number of stars should 
be: 
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n= 4352( 2 2) 
d\d 


ye eed pl 
12 \ 12 
= 922 stars 


The efficiency of the bursting charge when 
reacted should be: 


total moving energy of stars 
r= 


energy of the burst charge 


= J -mV?n [og 
2 y-l 

= +0.00153x(167.0)' «922 | HTx9-80 
(25-1 

= 0.0070 

=0.7% 


where y is the adiabatic thermal expansion coef- 
ficient of the burst charge gas. The gravitational 
acceleration g appears in the above equation 
because the force of explosives f was calculated 
on the basis of a force of 1 kg-weight x 1 meter, 
which is 9.80 joules. To convert (kg-weight x 
meters) into joules we must multiply by 9.80. 
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Chapter 9 — Discussion and Conclusion 


Tests of bursting warimono shells that were 
manufactured under various conditions were 
conducted on the ground, and the position of 
moving stars were analyzed by a photographic 
method. To make the analysis easy, ring star 
shells were employed as test shells. The results 
showed the quantitative relations among design 
factors, which were not clear until this study. 
These data were introduced into empirical equa- 
tions concerning the ballistics of stars, which are 
useful for designing the chrysanthemum shell. 


Convenient tables to give clear information 
concerning designing, characteristics of stars 
and shells were prepared. 


From the above effort the following infor- 
mation was obtained: 


1) The burning velocity of the stars is the same 
whether in the air or on the ground. There- 
fore, for designing the star, we can use the 
burning velocity on the ground. 


2) The characteristics of the star moving in air 
depend on the two values A, and b. The 
smaller the value of A, and the larger the 
value b, the better the characteristics of the 
star. The term A, is proportional to 1/6. 
Therefore, larger burning velocity and larger 
density 6 result in better star characteristics. 
The desirability of a large burning velocity 
has been commonly known, but the effect of 
larger density was first shown in this paper. 
And b is denoted as: 


Therefore, the larger the inside density, the 
better the star characteristics. If we design 
long burning stars, such as willow, it must 
be thought to be quite the opposite. 


3) From experiments, we find that air resis- 
tance to the motion of the star is propor- 
tional to the square of the velocity of the 
star. However, for very slow moving stars, 
future studies are needed. 
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4) On the bursting charge, the force of explo- 
sion f was obtained by calculation, and the 
relative value of the vivacity A was obtained. 
The relative values coincided closely with 
the reciprocals of the burning times of the 
bursting charge grains. Therefore, when 
planning the bursting charge, the value of A 
is determined by the burning time in air on 
the ground. On the calculation of the values 
of f it is better to compare with the values 
that will be obtained in the future by bomb 
tests. 


5) For the design, the value A should be as high 
as possible. Therefore, when using a burst- 
ing charge of low burning velocity, like po- 
tassium perchlorate powder, it is necessary 
to select a reasonable grain size. Specifically, 
slow burning powder may be useful with a 
proper grain size. However, it must be noted 
that too large a value of A causes ignition 
failure or breaking of stars. In this study, it 
is thought that the value of A should be from 
3 to 5. 


6) Many layers of pasting paper increase the 
efficiency of the bursting charge. However, 
if there are too many layers, it has the re- 
verse effect. This effect may be caused by 
unnecessary consumption of energy to break 
the thick-pasted shell. 


7) The ratio of the break strength of the Japa- 
nese paper and that of the Kraft, which was 
tested by a tensile strength instrument, was 
slightly different from that from the bursting 
experiment. However, when there are no 
values from bursting experiments, the tensile 
strength test will be effective in deciding the 
value of a. 


8) In the above experiments the differences be- 
tween ring and full star shells were not clear 
because there were too few full star shells. 
This point may be clarified in the future. 


9) In this study the loading density was set as a 
constant. This problem remains for future 
investigation. 
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10) The use potassium perchlorate bursting 


charge in place of potassium chlorate burst- 
ing charge, which has been though to be 
much more dangerous than the former, was 
clarified. 


11) An example of designing a chrysanthemum 


shell is shown. It shows that the fundamen- 
tal method to plan the chrysanthemum 
shell, which has previously been made only 
by application of some traditional ideas, has 
been developed to use physical or chemical 
science. 
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